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ABSTRACT 



Aims. A thorough study of the effects of mass loss on internal and surface abundances of A and F stars is carried out in order to 
constrain mass loss rates for these stars, as well as further elucidate some of the processes which compete with atomic diffusion. 
Methods. Self-consistent stellar evolution models of 1.3 to 2.5 Mq stars including atomic diffusion and radiative accelerations for 
all species within the OPAL opacity database were computed with mass loss and compared to observations as well as previous 
calculations with turbulent mixing. 

Results. Models with unseparated mass loss rates between 5 x lO"''' and 10"'^ Mgyr"' reproduce observations for many cluster AmFm 
stars as well as Sirius A and o Leonis. These models also explain cool Fm stars, but not the Hyades lithium gap. Like turbulent mixing, 
these mass loss rates reduce surface abundance anomalies; however, their effects are very different with respect to internal abundances. 
For most of the main-sequence lifetime of an A or F star, surface abundances in the presence of such mass loss depend on separation 
which takes place between log AM/M, = -6 and -5. 

Conclusions. The current observational constraints do not allow us to conclude that mass loss is to be preferred over turbulent mixing 
(induced by rotation or otherwise) in order to explain the AmFm phenomenon. Internal concentration variations which could be 
detectable through asteroseismic tests should provide further information. If atomic diffusion coupled with mass loss are to explain 
the Hyades Li gap, the wind would need to be separated. 
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1 . Astrophysical context 



As instruments become more sophisticated and precise observa- 
tions are readily made available for a growing number of chemi- 
cal elements, additional constraints are steering the evolution of 
stellar models. The inclusion of atomic diffusion and radiative 
accelerations into the standard stellar evolution model resulted in 
some early success in descr ibing abundance patterns for AmFm 
stars dTurcotte et aU figgSaV In these models, particle transport 
within the radiative zone was calculated from all physics known 
through first principles. However, the calculated surface anoma- 
lies were greater than the observed anomalies, and it was deter- 
mined that additional transport processes were flattening the sur- 
face abundance patterns. With the addition of turbulent mixing , 
the Montreal g roup's stellar evolutio n code ([Richer et al. , 200C ; 
Richard et al 1 120011: iMichaud et all l2004 120051: iTalon et al 
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2006h was able to explain and reproduce the particular abun 



dances of many cluster and field stars with a single tunable pa- 
rameter: the mixed mass. Nonetheless, other studies have sug- 
gested that mass lo ss could also reduce the predicted anomalies 
to observed levels jMichaud et all 119831: iMichaud & Charlandl 



ll986l:lAlecianL ll9m Since these studies considered static stel- 
lar models incorporating a limited number of species, a rigorous 
investigation of the effects of m ass loss on stellar chemical evo- 
lution is warranted. Previously, IVauclair & Charbonnell (Il995h 



have introduced mass loss to reduce the effect of Li settling in 
Pop II evolutionary models. 

AmFm stars (7000< T^ff < l OOOOK) are slowly rotating 
(v,ot < lOOkms \ lAbtl l2000h . non-magnetic stars of the 
main-sequence (MS). They are interesting candidates for test- 
ing evolutionary models because they lie within the temper- 
ature range for which the depth of the surface convection 
zone varies quite rapidly. Chemical separation within the sta- 
ble radiative zones of these stars generates surface abundance 
anomalies within timescales that depend strongly on the depth 
of the surface convection zone, although the exact depth at 
which the separation occurs is still debated. The original ex- 
planation for these chemicaUy peculiar stars stipulated that 
the separation occurred immediately b e low the surf ace H con- 
vection zone dWatsonL Il97ll: lAlecianL Il986l Il996l) : however, 
with the inc l usion of turbulent mixing, more recent models 
(iRicher etalil2000[) have suggested that separation occurs much 
deeper in the star (r > 200 000 K). Though the second scenario 
has had success not only with AmFm stars, but also with Pop II 
(Ri chard et al .. 2005; Korn et al., 2006) and Horizontal Branch 
(HB) stars JMichaud et al.> ,2007. i2008i) . it is still premature 
to accept turbulence as the sought-after macrospcopic process 
since all models with turbulence necessarily involve at least one 
adjustable parameter Turbulence models often implicitely as- 
sume a link between turbulence and rotation; however, even the 
most slowly rotating AmFm stars have anomalies which are sig- 
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nificantly smaller than those obtained with atomic diffusion only 
models. This suggests the presence of a competing process even 
in non rotating stars. 

Typical anomalies on the surface of AmFm stars include 
overabundances of iron peak e lements, as well as underabun- 
dances of Ca and/or Sc (see ICavrel et alj 1199 ll for a more 
complete description). Recent studies have obtained abundance 
determinations for numerous A and F stars of open clus- 
ters fOT_a_jiumberofche ( Burkhart & Coupry, 
200tf, 'Hui-Bon-Hoa', '200(/, 'Moniei*, '2005; Fos sati et al.. ,20071: 
Gebran et al., 2008; Gebran & Monier, 2008) in an effort to con- 
front relatively well constrained stars to current evolutionary 
models. The advantage of observing cluster stars is that they 
generally have the same age and initial metallicity, which greatly 
facilitates a comparison with models. 

At the cool end of the Fm star domain, the well doc- 
ume nted Li dip, first ob s erved in the Hyades open clus- 
ter (iBoesgaard & Tripiccoi 11986'), has challenged theoretical 
astronomers for decade s (Michaud, 1986; Vauclair, 1988; 
iMichaud & Charbonneaui ri99k iTalon & CharbonnelL i2005l) . 
The lithium abundance has also been obse rved in many other 
open clusters (e.g. [B alachandran 1995; iBurkhart & Couprvl 
l2000t |Anthony-Twaro2 et al. 2009). The Be abundances in 
these stars (BoesHaard& King, 2002^ iBoesgaard et al.L 120041; 
iRandich et all l2007l) provide additional constraints on parti- 
cle transport. Recently. iTalon & Charbonne l (2003, 2005) have 
modeled shear turbulence induced by differential rotation and 
mixing induced by internal gravity waves in order to describe 
both the hot and cold side of the dip. Other models have also 
explored the effects of horizontal /i-gradients on rotationally in- 
duced mixing in these stars (I'alacios el al. 2003). However, the 
potential effects of atomic diffusion, more specifically of radia- 
tive accelerations, in competition with mass loss, were not fully 
investigated. 

At some level, mass loss is present in all stars, and 
it is important to quantify its effects on observed abun- 
dances. Unfortunately, for A and F stars in particular, the 
mass loss rates are not known. In O and B stars, the ra- 
diatively driven winds produce mass loss rates as important 
as lO^'^M^vr ' dLamers & CassinelUl Il999l) . In colder stars of 
types G and K, winds are driven by active coronas. The best 
known example is the solar wind whic h has a mass loss rate of 
2 X lO '-^Mayr-i ( Feldman etal]ll977h . For intermediate stars, 
our understanding is at best nebulous. IXbbott (1982) suggests 
that radiative accelerations are too small in stars with reff< 
lOOOOK for radiatively driven mass loss to be significant. On 
the other hand, the thinning of the surface convection zone in 
F and particularly i n A stars mi ght be too important for solar 
type winds to exist (IParker||l960). Are winds of A and F stars 
driven by radiation, coronal heating, both or neither? What are 
the expected mass loss rates? Both of these questions remain 
unanswered, and answering one could shed light on the other 

Observ ational constraints on A and F star mass loss rates 
are limited. l Lanz&Catalal([T99l as well as lBrown et al. 
gave an upper hmit of lO"'*'M0yr 'for main seque nce A stars. 
Asy mmetries in Mgll spectral lines of Sirius A led iBertin et al.l 
(11995) to conclude that mass loss is present with a rate be- 
tween 10 '^Moyr ' and 10 '^Moyr On the theore tical side, 
the radiatively driven wind model of iBabell (Il995h suggests 
a mass loss rate of 10~'^Moyr ' for A stars. However, ac- 
cording to his results, only heavier elements are evacuated by 
the radiative field. S imilarly, Michaud et al. (1983, see also 
IMichaud & Charlandl [l986) suggested that mass loss rates be- 
tween lO^'^^Moyr"' and lO^'^Moyr"' could satisfy observa- 
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Fig. 1. H-R diagram for all the models shown in Fig.|4] Though 
all models were calculated from the PMS to the bottom of the 
subgiant branch, the complete tracks are only shown for the 1 .7 
and 2.3 Mq models. All other models are shown only on the 
main-sequence, which, for the purpose of this plot, is consid- 
ered to span from the time at which 1 % of core H is burned to 
the time at which 95% is burned. The dotted line represents evo- 
lution on the PMS until diffusion starts (i.e. until a radiative zone 
appears). 



tional constraints from CP star surface anomalies. Given the 
large disparity in values, the mass loss rates used in this study 
will be constrained stricly by surface abundance variations (see 
also Sect. 4.2). 

In this paper, we consider mass loss in non rotating stars. In 
these stars, mass loss is arguably the only macroscopic process 
competing with atomic diffusion within the radiative zones. We 
will start with a brief description of our stellar evolution code in 
Sect.|2]after which we will discuss the method for calculating ra- 
diative accelerations for lithium, beryllium and boron (Sect.O. 
In Sect.|4]we will describe the treatment of mass loss. In Sect. 5 
we will discuss its effect on internal structure and surface abun- 
dances as the models move along the main-sequence and the 
subgiant branch. In Sects. 6 and 7 we will compare our models 
to turbulence models and observations respectively. In Sect.[8] a 
brief overview of the main results will be followed by a discus- 
sion on how asteroseismology could help decipher the effects of 
advective processes including meridional circulation from those 
incurred through turbulent processes caused by differential rota- 
tion. 

2. Calculations 

This pap er is part of the Montreal stellar mode l development 
project (iRichard et alJ200TllTurcotte et al.ll998b and references 
therein). The models were evolved from the initially chemi- 
cally homogeneous pre-main-sequence (see Fig.[T]) with the so- 
lar abundance mix listed in Table 1 of Turcotte et al. (1998^. 
The transport of chemical elements in 1-D (one di mension) 
is solved within the basic framework established by iBurgersI 
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(Il969l) . The chemical transport equation considers all effects 
of radiative accelerations and atomic diffusion, and it is solved 
for 28 chemical elements and isotopes at every mesh point (the 
number of points varies from 800 to about 2800) for each time 
step. Atomic diffusion is allowed to operate as soon as a ra- 
diative zone appears. The radiative acelerations (giad) and the 
Rosseland mean opacity are continuously updated; the treatment 
of chemical transport is thereby fully self-consis tent. The ato mic 
diff usion coefficients were ta ken from Paquette et al.l ( Il986l see 
alsolMichaud & Proffit jl993l) . The Krishna S wamy T-t relation 
dKrishna-Swamv, 196^ was imposed as the boundary condi- 
tion in the atmosp here (this c hoice was motivated by results pre- 
sented in^VandenBerg et al.'2008). Semiconvection was included 
as described in Richard et al. (2001), following iKatol (Il966h . 
ILanger et al.1 ( ll985i) . and .Maedena997i) . 

The chosen values of the mixing length parameter and initial 
He fraction are respectively a = 2.096 and Yq - 0.27769 (see 
Model H of lTurcotte et al.ll998bl) . which are calibrated by fitting 
the current solar radius and luminosity. We chose Zq = 0.01999 
as the initial mass fraction of metal^ . Some models were also 
calculated for Zq = 0.01 and Zq = 0.03. 

These are the first fully self-consistent stellar models which 
include mass loss. Models were calculated from 1.30 
to 2.50 Mq. The mass loss rates considered vary from 
1 X lO""'M0yr"' to 1 x lO^'^Mgyr"'. Our treatment of mass 
loss will be further discussed in Sect.l4.1l 



3. Radiative accelerations 

Rosseland opacities as well as radiative accelerations are contin- 
uously computed for all 28 species as the relative concentration 
of each species varies with time. For all elements included in 
the OPAL database felesias & Rogersl 1 19961) . the radiative ac- 
celeration calculations are carried out using direct summations 
over the actual spectrum (i.e. opacity sampling, see lRicher et al.l 
11998 ). At large optical depths (where the diffusion approxima- 
tion is valid), the radiative acceleration, grad(^), of an element A 
at a radius r in a star may be approximated by: 



5rad(A) 



1 Lr KR 

4nr^ c Xa 



Jo Kuitotal) 



'P(u)du, 



(1) 



where 'P(m), the normalized black body flux distribution, is given 
by: 



'P{u) = 



15 



47r4(e«- 1)2' 



(2) 



The radiative luminosity at a radius r is U^^^, Xa is the mass 
fraction of A, kr is the Rosseland opacity, /<-„(total) and Ki,{A) are 
respectively the total opacity and the contribution of A to the 
total opacity at the frequency u defined by: 



u = hv/kT. 



(3) 



Since the competition for photons between element A and all 
other elements present in the plasma determines the value of 



' lAspIund et all ( l200l l2009h have proposed a downward revision 
of the solar abundances of some metals; however we have chosen to 
keep the previous abundances until t heir determ inations are reconciled 
with helioseismology (Delahaye & Pinsonneaul t 2006; Basu et al. 2007 
and references therein). The abun dance of O, the third most abundant 
eleme nt, is particularly uncertain dCaffau et all 120081 : iDelahave et al.L 

Eoia). 



giad{A), atomic data is required for all species in order to com- 
pute grad for any given element. 

The corrections f or re distribution of mom entum are from 
iGonzalez et aH (1 19951) and lLeBlanc etaP (l2000h . 

3.1. Radiative accelerations for Li, Be and B 
3.1.1. Methods 

Li, Be, and B are not included in OPAL since they are not 
significant in the calculation of Rosseland mean opacitjQ. We 
have nonetheless computed the radiative accelerations in a man- 
ner which is consistent with OPAL spectra. The various cor- 
rections determined by IRicher et al.l d 19971) are al so included. 
The at omic data required for LiBeB are taken from lWiese et al.l 
(119661) . Since these elements are not sufficiently abundant for 
pressure broadening to be important (the lines are never satu- 
rated), only oscillator strengths are required to compute opaci- 
ties. 

The calculation of grad for these elements is delicate because 
of two factors which can lead to large fluctuations: (1) non op- 
timal frequency sampling and (2) random background changes. 
Both problems arise from the fact that LiBeB not only have very 
few lines contributing to their giad, but these lines are also very 
narrow. 

(1) The fluctuations due to sampling not being sufficiently 
refined has an effect on Li as well as on the background (see 
Fig.|2]i. The Doppler width of a ^Li line is given by: 



Am 



2kT 



0.5 



5.1 X lo^^r"-^ 



(4) 



where u is the adimensional energy difference between the upper 
and lower levels of the transition, M is the mass of the element, T 
is the local temperature and T5 - T/IO^ K. In this case, one typ- 
ically has M ^ 5 around the most important temperature range, 
log T ^ 5.3, so that 



AMDopplei- - Am^ 



^7x10" 



(5) 



which is to be compared to the frequency grid interval, Amg = 
20/10'* = 2x10"^, of the OPAL spectra. Using opacity sampling 
for Li would imply that the line center would sometimes occur 
within Am/5 of a grid point and sometimes up to 15Amo away. If 
the sampling point occurs in the far tail of an important Li line, 
the §i-ad(Li) value would be much smaller than if it occurs in its 
center. 

Likewise, the sampling is not refined enough to reproduce 
all features of the background spectra (see Fig. 3 of Riche r et"al1 
1998). Variations can therefore result from overlooking an im- 
portant background line which occurs near a lithium line. 

(2) The frequency sampling grid is a function of u and not 
of V. Therefore, when one considers an adjacent grid point of 
different T, the background, as a function of u, changes, and 
the consequent random variation in background affects the flux 
available for a specific lithium line. This is largely due to the 
narrowness of most lines. For example, in the inset of Fig.|2] the 
Li line does not encounter any important Fe line for that specific 



2 In lTurcotteetai] j 1998b) and Turcott e et al.l ( Il998al) . grad(Li) and 
ftad(Be) were app roximated using formulae derived previously by 
Richer & Michaud ( 1993). These are however less accurate than giad 
calculated using direct summations over the spectrum throughout stel- 
lar evolution (using Eq. [T]), which are descr ibed here and were used in 
iRicher et alj ( I2OOQ) and iRichard et alJ (1200 ll) . 
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Fig. 2. Opacity spectra for Li, Fe as well as the total opacity in cm^g"' at a depth where log T = 5.3 and log p = -4.6. The u interval 
covers the range through which most of the flux passes. The spectrum is dominated by Fe lines. The inset on the right shows a zoom 
of the area which hosts the Li line which contributes the most to its g^d- 



i p,T) t able; however, in the inset of Fig. 1 of iRicher & MichaudI 
2005h . which shows the same u interval for a different (p,T) ta- 
ble, the Li line is overlapped by a strong Fe line, and the available 
flux is consequently reduced. 

One can alleviate these problems by combining a modified 
version of the opacity sampling method with the known position 
of all the line centers for LiBeB. In order to reduce the fluctu- 
ations, while preparing the spectra of, say, Li, for the calcula- 
tions of gra(j(Li), it was assumed that each Li line was spread 
uniformly over the Auq interval in which the Li line center oc- 
curs (i.e. the Li lines become square functions of width Auq and 
so it becomes impossible for the sampling to overlook them). 
This will leave fluctuations caused by random variations of the 
opacity background when the line center moves from one Amg 
interval to another (2), but the variations due to non optimal sam- 
pling (1) will be reduced significantly. 

Note that if the frequency sampling is refined (e.g. 10^ as 
in OP data) in order to better represent the background spectra, 
errors due to the inexact position (in frequency) of each back- 
ground feature remain. The line center positions for Fe, the main 
contributor to the total opacity in this region, are only known to 
about 1%. 

The resulting error bars on gi-^d are discussed in Sect. 3 of 
iRicher & MichaudI (1200 5). As an example, for population I stars, 
there is a factor of 2 uncertainty for grad(Li)- In Pop II stars, the 
uncertainty is much less important (e.g. only 3% for a star with 
Z = 0.0001). 

The same errors and limitations should be expected for scan- 
dium around the minima of its ^lad since the radiative accelera- 
tions in these regions are computed with only a few narrow lines. 
Unfortunately, Sc is not included in OPAL d ata, and its grad must 
be ca lculated through alternate methods (iLeblanc & AlecianI 
I2008h . 

3.1.2. Results 

The radiative accelerations for Li and Be are shown in Fig.[3]for 
a star of L55 Mq with a mass loss rate of 1 x lO"'"'M0yr"' as 
well as for a star of 1 .47 Mq with no mass loss. For both models, 
the radiative accelerations for Be are below gravity throughout 





-1 
^ -2 

a 

M) 
O 

" -1 

-2 





' 1 ' 




' 1 ' 


Be : 


I — \ — \ — 


H h- 


H 1 1 1 1 1 


1 1 h 


1 r^: 

Li 




, 1 , 











5.0 


5 5 
log t (K) 


6.0 





Fig. 3. Variation of radiative accelerations with temperature for 
Li and Be at lOOMyr (dotted fine), 700 Myr (dashed fine) and 
L3 Gyr (long dashed line) in a L55 model with a mass loss 
rate of 1 x lO^'-'Moyr"' as well as at 100 Myr (soHd line) for 
a L47 Mq model with no mass loss. The curves extend to the 
bottom of the surface convection zone, on the left. Acceleration 
is normalized to local gravity. The horizontal arrows indicate the 
total movement of the bottom of the surface convection zone 
(BSCZ) for models of 1.40, 1.43, 1.45, 1.47 and 1.55 Mq from 
top to bottom (see text below). 



the radiative zone, and throughout evolution. However, §rad(Li) 
has a peak that approximately reaches gravity below the surface 
convection zone for the 1.55 model. Up to log T ^5.4, grad 
is seen to be almost solely temperature dependant as the curves 
for both models overlap. 

In the 1 .47 model, g,ad(Li) has about the same value as 
in die 1.55 Mq model for log T > 5.4, but, for log T < 5.4 it is 
reduced by the competition of Fe, which absorbs most of the flux 
where the Li lines are important. This occurs when mass loss is 
small enough for Fe to accumulate in and above this region, as 
is the case for the 1.47 model (or in the 1.50 model with 
a mass loss rate of 10 '"'Moyr"' shown in Fig.|5]l. On the other 
hand, grad(Be) does not vary with Fe abundance likely because 
most Be lines, particularly the ones which contribute the most to 
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its gmd, lie outside of the frequency interval which is dominated 
by Fe lines. 

The horizontal arrows in Fig. |3] show the total movement of 
the bottom of the surface convection zone (BSCZ) for 4 mod- 
els without mass loss (1.40, 1.43, 1.45 and 1.47 Mq), as well as 
for a model of 1.55 Mq with a mass loss rate of lO^'^'Moyr"'. 
The interval for which the BSCZ is plotted spans from 50 Myr 
to 625 Myr for all models except for the 1 .47 Mq, model which 
goes from 50 Myr to 425 Myr (the simulation's last converged 
model). For most models, the depth of the BSCZ does not vary 
much between 50 and 625 Myr (most of the main-sequence life- 
time before the age of the Hyades cluster); however, for the 
1 .47 Mq model, the BSCZ moves significantly and extends over 
an interval in which giad(Li) also varies significantly. 

If we compare the curves for the 1.55 M© model with 
those shown in Fig. 1 of Richer & Michaud (1993), we first 
note the similarity in the temperature dependence of the 
curves. For Li, we have verified that they are nearly equal for 
logT > 5.4. The maxima occur at very nearly the same tem- 
perature as well. However, the maximum values of grad(Li) 
and ^rad(Be) are respect ively about 2 and 6 times larger in 
[Richer & MichaudI (Il993h . Furthermore, their results show very 
smooth curves compared to the many variations seen in our cal- 
culations. 

Although there is a slight difference in stellar mass ( a 
1.54 Mq model is shown in Fig. 1 of Richer & Michaudlll993h . 
this cannot account for the relatively large differences close to 
the maximum. There are nonetheless a few other explanations. 
Our calculations were carried out with integrations over com- 
plete OPAL spectra for all species at each time step and each 
mesh point. Since these spectra were not av ailable in 1993, 
the calculations of [Richer & MichaudI (Il993h were done us- 
ing averaged spectra whose frequency dependence (Eq. [18] of 
[Borsenberger et al. 1979) did not include the frequency depen- 
dence of Fe lines, particularly near T ^ 2 x 10^ K, where grad(Li) 
is affected the most. From Fig.|2] one can see that Fe lines dom- 
inate the region where Li lines are strongest. As Fe lines which 
occur near Li lines were absent in the old opacities, the available 
flux for Li lines is very different (see Sect. 13.1. Ti l. 

Moreover, the evolution of individual metal abundances and 
its impact on local opacity are not included in the calculations 
of [Richer & M ichaud (1993). This is particularly important in 
AmFm stars since heavier metals such as iron and nickel tend to 
accumulate below the surface convection zone when mass loss 
is not too strong. 



locity due to the wind: 



v„(r) 



M 



Mt 



(6) 



where p is the local density, m,- is the mass interior to r and 
M < 0. It should be pointed out that the m^/M, factor in this 
equation has often been overlooked, but it will be shown to be 
a consequence of mass conservation (see Sect. 14.1. Ti l. Even for 
small mass loss rates, the p"' dependence leads to large veloci- 
ties in the outer layers of the star. To avoid numerical instabilities 
due to the cancellation of two large quantities in the convective 
envelope (wind velocities and large turbulent velocities used to 
enforce convective mixing) and in order to have a Neumann sur- 
face boundary condition (no flux), mass loss was implemented 
as described in ICharbonneaul (119931) . The main physical con- 
siderations are that the surface convection zone (SCZ) be fully 
mixed, and that the atmosphere be mixed with it by overshoot- 
ing. Therefore, the ejected mass has the same composition as 
the SCZ. The mechanism by which mass is ejected from the star 
becomes irrelevant so long as the correct amount of mass is re- 
moved from the star For each chemical species, the resulting 
transport equation is: 

P^ = -V-[-pDVlnc + p(U + UJc] 

+p(S „uc + S ,,)c, (7) 

with a Neumann condition imposed at the surface and with U„ 
and Sk defined as: 



VvvCr under the SCZ, 
in the SCZ; 



under the SCZ, 
= ^ £intheSCZ. 



(8) 



(9) 



Here, c is the time and depth dependent composition, D the total 
diffusion coefficient, U the advective part of the atomic diffusion 
velocity, the interior velocity due to the wind, Mqz the mass 
of the SCZ, M the mass loss rate, S„iic a source/destruction term 
linked to nuclear reactions and S „ is a sink term linked to mass 
loss. The mass is simply removed from the surface convection 
zone which, given the assumed mixing, is equivalent to losing it 
through the stellar surface. The mass of the star is also continu- 
ously updated so that all quantities that depend on stellar mass 
are correctly calculated. 



4. Mass loss 

We first discuss how to include mass loss in evolutionary models 
(Sect. 14. it . In Sect. 4.2, we discuss a few theoretical analyses of 
coronal and radiatively driven winds followed by a brief look 
into separated win d£| (Sect. 4.3). 

4. 1 . Treatment of mass loss 

Mass loss is assumed spherical and chemically unseparated. If 
we simply apply mass conservation arguments, the net result of 
mass loss is the appearance of an outward flowing interior ve- 



^ In this paper, mass loss is separated when the abundances in the 
wind are not the same as in the photosphere. On the contratry, mass loss 
is unseparated when the abundances in the wind are the same as in the 
photosphere. 



4.1 .1 . Mass flux equation 

Mass loss has the effect of ejecting (or peeling) the outermost 
layers of a stellar model. The model must then be reconverged 
with a slightly reduced mass. This scenario can be described 
within the formalism of operator splitting. If a mass AMioss is re- 
moved in a given time step Af, one needs to reconverge a model 
of mass Ml = Mq - AMioss at time f i = fo -H Af with the composi- 
tion of the star kept unchanged as a function ofnir, the indepen- 
dent variable which is defined as the sphere of constant mass, mr. 
Mass loss does not modify the abundance profiles with respect 
to nir, although there is a change with respect to the stellar sur- 
face. For example, an abundance peak at the point Mq - AMioss 
(or AMioss below the surface before peeling), finds itself at Mi 
(or at the surface) after the peeling but the mass interior to this 
point, m,., has not changed. 

To implement mass loss in a stellar evolution code, one can 
imagine the solution process to be broken down into 2 steps. In 
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a first step, the model can be converged at a time t[ as if there 
were no mass loss, with the composition changes due to nuclear 
reactions and diff'usion processes. Then, the grid at fi is reinter- 
preted as corresponding to a mass Mi = Mq - AMioss with all 
variables expressed as a functiorQof nir/Mt except for the com- 
position which is kept unchanged as a function of m,-. The star is 
then reconverged a second time at f i with a mass Mi and with the 
concentration profiles which include the effect of mass loss. As 
more mass is lost, the model inches toward less massive struc- 
tures, and so evolution progressively corresponds to that of a 
lower mass star. This is a correct description of mass loss which 
in practice is as accurate as the operator splitting procedure is 
accurate. 

However, going from the nir/Mo to the nir/Mt, grid while 
keeping c, as a function of nir implies interpolating. In practice, 
performing interpolations on c, during computations can intro- 
duce unwanted numerical diffusivity. It is relatively easy to show 
that these interpolations on c, resulting from mass loss can be 
avoided by introducing a local mass flux into the conservation 
equation, while keeping the grid constant (i.e. the grid points 
have the same values of m,./M, in spite of M, varying from Mq 
to Ml). 

Let us consider the concentration on a grid point (mr/M,) at 
ti as a function of the concentration on the same grid point at to. 
The grid points shift on the m, axis. Then 

dci 

Ci (nii-i) = Ci (m,i)) + - — Anir (10) 
onir 

which is a simple Taylor series development and where 

Anir — m,-\ — m,-Q. (11) 



By definition, m,-i is the position on the axis, of the grid point 
^ after a AMioss mass loss so that one may write: 

niri niii) + Am,- niro 



Ml Mo - AMioss Mo 
which can trivially be solved to give: 

^"^1- " -Tr(AMioss). 
Mo 



(12) 



(13) 



Now we substitute equation ( [13] ) in equation ( fTOl l. replace AMioss 
by -M X At and simplify: 



, , , , dci dM, niro . ^ 
drrir df Mo 



(14) 



The second term on the right has the same effect on c, as the 
introduction of a flux terrrQ (caused by mass loss Xj^) in a con- 
servation equation. The flux Anpr^Uw given by equation 



- 4-nprU», = 



, dM, m,- 

,ljT » I 



dr M. 



(15) 



must then be introduced in the conservation equation in order to 
take into account the effect of the nir shift caused by the peeling 
of surface layers while keeping the same grid as a function of 
nir I when M, changes due to mass loss. This simple argu- 
ment justifies the introduction of a mass flux to model the effect 



Here Mj = M, at ti and similarly Mq = M, at fo- 
^ Often abusively called a "wind" term in this context. It seems 
preferable to restrict the use of "wind" to the region above the pho- 
tosphere. 



of mass loss on element separation in stellar atmospheres or e n- 
velopes (I Vauclaiii [1971 iMichaud etaLl 1 1 9831: lAlecianL 1 1 996l) . 

For smaller mass loss rates, structural effects of mass loss are 
often negligible to the extent that only the shift of c, remains and 
mass loss may be viewed as a mass flux going through a star of 
constant mass (i.e. it is not necessary to change the stellar mass 
during evolution calculations). Consider the case of the present 
Sun. Assuming that it is constant in time, its current mass loss 
rate of 2x lO^'^^Mo/yr leads to a loss of 10""* Mq up to the Sun's 
age which for most purposes corresponds to negligible structural 
changes. This is only ~ 1/300 of the mass of the superficial 
convection zone. 

Since we are not certain about the nature of the winds at the 
surface of A and F stars, the present models do not take into 
account any energy dissipation which would be required to pro- 
duce these winds. 



4.2. Stellar winds of A and F stars 

This study will use observed surface abundances to constrain 
mass loss since the stellar winds associated with A and F 
stars are not well known. The winds could be radiative, coro- 
nal, a combination of both or even completely negligeable. 
Comparisons with stellar wind models are difficult. Indeed, even 
the hottest A stars maintain a thin surface H co nvection zone, 
and stars as early as A7 (and pos sibly earlier, see iNeff & SimonI 
2008; Si mon & Landsmanll997h can support active coronas and 
chromospheres which could harbor solar type winds. It is also 
plausible that both mechanisms act simultaneously. A few prop- 
erties of coronal and radiatively driven winds for A and F stars 
will be reviewed in the two following subsections, in so far as 
they relate to chemical separation. 



4.2.1 . Coronal winds 

With simple physical considerations, it is possible to obtain an 
approximate value of the mass loss rate above which coronal 
winds are necessarily unseparated. For a spherically symmetrical 
mass loss, the wind velocity (Eq. ||6l with Mr = M,) may be 
compared to the maximal downward diffusion velocity (given by 
the gravitational settling velocity without any contribution from 
radiative accelerations). Equating the two gives an evaluation of 
the maximum mass loss that allows for separation to occur One 
may then write: 



-dM, . 

Anr^p kT 



which can be rewritten as: 



-1 dM, 2.4x10 



-15 



M, df 



AiT\-' 1 



(16) 



(17) 



{T^ - T/IQ^ K) where A, and Z,- are the atomic mass and atomic 
number respectively of element / and where the Coulomb term 
in the calculation of Djp is replaced by an approximate value. 
This applies both in the atmosphere and in the outer parts of 
the wind solution. For a solar type wind, assuming an isother- 
mal corona of T = 10^ K, one obtains a limiting mass loss rate 
^ lO"'^M0yr"'(see Michaud et al. 1987). The dependence on 
T is valid when the gas is fully ionized. In the cooler atmo- 
spheres of early A stars (say at lO'* K), He is mainly neutral, nev- 
ertheless the limiting mass loss rate is approximately the same 
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in the atmosphere as in the corona since the diffusion coeffi- 
cients are larger by a factor of 200 - 300 than those of ionized 
He (Michaud et al. 1978). The separation observed in the solar 
wind (M ^ lO^^^Moyr-i) agrees roughly with this analysis. As 
previously mentioned, observations show that elements which 
are neutral at T - 10"^ K are ^ 2 - 3 times less abun dant in 
the corona than in the photosphere ('Meyejl 1 19851 Il996l) . Most 
agree that the separation occurs somewhere in the solar chromo- 
sphere immediately above the photosphere and could therefore 
have ar i impa ct on obs erved p hotospheric abundances. For in- 
stance, iGeiss & Buergil (Il986i in their model 4) find that matter 
which arrives at 1 AU has the same He to H ratio as at the wind's 
base. Since there is no net separation within the solar wind, it is 
precarious to infer that all coronal type winds are chemically dif- 
ferentiated. Moreover, wind structures in A star coronas would 
probably be quite different since many metals that are neutral in 
G and F star photospheres are ionized in A star photospheres. It 
is thus very difficult to compare the solar wind to other coronal 
winds, particularly for hotter stars. 



4.2.2. Radiatively driven winds 

iBabell (1 1995b found that for stars within 8000 < T^s < 
14000 K, all radiatively driven winds must be fully sepa- 
rated. In these stars the Coulomb coupling is not sufficient 
to redistribute the momentum acquired by the heavier, radia- 
tively accelerated elements onto t he bound, more abundant H 
and He (iSpringmann & Pauldrachi [T992: Owocki & Puis. 2002; 
iKrticka et al.L I2003h . Since only metals are ejected from the 
star, the mass loss rates are much smaller: between 10 '^M© 
and 10"'^ Mq. Ho wever, the m ulticompon ent hy drodynamical 
model put forth in lBabell (Il995h and lBabell d 19961) only consid- 
ers an average metal, rather than solving for each metal indi- 
vidually, therefore metal-specific mass loss rates are not known. 
Furthermore, the interaction of radiatively driven winds with 
magnetic fields as well as convection is still poorly understood. 
This is particularly important for A stars, and the consequent un- 
certainties require us to be cautious before constraining our anal- 
ysis with these results. For cooler stars, such as F stars, radiative 
accelerations are not believed to be able to generate significant 
mass loss (Abbott, 1982). 

lUnglaubl (l2008h found that radiatively driven winds of sdB 
stars must be separated (accelerated metals cannot drag H and 
He) if the mass loss is smaller than lO^'^Moyr"' (10 times larger 
than the approximate value obtained in Sect. l4.2.TT i. However, al- 
though the author's calculations are quite thorough, they are not 
complete in so far as some other poorly understood factors could 
play a significant role in determining wind properties. How does 
convection or magnetic fields, particularly flux tubes, affect the 
wind structure and velocity? Furthermore, the omission of line 
shadowing in the calculations could have a significant impact on 
the author's result^ In fact, the author stipulates in Sect. 6.2 of 
his article that including line shadowing could diminish gi-^^ by 
a factor of 100 for the stronger photoshperic lines, thus leading 
to an overestimation of M by a factor of 10. 



* Line shadowing occurs when the wind velocity is not sufficient to 
Doppler shift the wind's line centers away from flux attenuated photo- 
spheric line centers, thus reducing radiative accelerations. This is often 
true until the wind reaches the sonic point. 



4.3. Unseparated vs. separated mass loss 

The object of this study is to constrain the effects of mass loss 
solely via observed abundance anomalies. To do so, we use sim- 
ple wind models in order to minimize the arbitrariness of the 
analysis. Accordingly, most calculations were done assuming 
simple unseparated winds, although a few calculations were also 
carried out assuming separation in the wind in order to assess po- 
tential effects. Three cases of separated winds were considered: 
(1) only metals are ejected, (2) the separation mimics the solar 
wind with H treated as a high-FlP (First Ionization Potential; 
lMeveiilT985h element and (3) the separation again mimics the 
solar wind but with H as a low-FlP element. 

In case 1, all metals are ejected with the same composition as 
the stellar surface, while H and He remain bound. For this sce- 
nario, mass loss rates were varied from 10"'^ to 5 x 10"'^ M^yr^^ 
in order to account for the fact that only metals are leaving the 
star (around 2% of the superficial mass fraction). 

The other scenarios (2 and 3) consider c hemic al separation 
in the solar wind as established by M eve^ d 19851) . who found 
that elements with a FIP smaller than 9eV were approximately 
4 times more abundant relative to hydrogen in the corona than 
in the photosphere, while higher-FIP elements, including hydro- 
gen, kept their photospheric abundances. Although this scenario 
is generally favored. M ever. (1996) questioned his own results 
a decade later by implying that instead of having overabundant 
low-FIP elements in the wind. higher-FIP elements, including H. 
could be depleted in the coroneQ. Both configurations are investi- 
gated: case 2 has H as a bound high-FIP element, and case 3 has 
it as a low-FIP element. Our approach was to divide all elements 
into two groups: low-FIP elements (below 1 1 eV) and high-FIP 
elements (He. C. N. O. Ne. CI and Ar). All low-FIP elements 
were depleted 4 times faster than high-FIP elements 

Numerically, in cases (2) and (3). the destruction term in 
Eq. Q was multiplied by a factor of 4 in the SCZ for all low- 
FIP elements. A weight term, which was continuously updated 
as concentrations changed in the SCZ. was added to the denom- 
inator for normally depleted high-FIP elements to account for 
the fact that their destruction is not. in this case, proportional to 
the total mass loss rate multiplied by individual concentration 
(Eqs. IQ and ||9)) since their relative concentration in the photo- 
sphere is not the same as in the wind. The interior wind velocities 
are not affected since the v„ind term does not depend on relative 
concentrations, but simply on the mass loss rate. 

5. Evolutionary models 

In Fig.m the evolution of T^s, L, \ogg, the depth of the surface 
convection zone as well as the surface Fe abundance are shown 
for a number of models. The chosen masses were selected to 
span the observed T^tf range associated wit h AmFm stars and Li 
gap stars (lPrestonll974tlBoesgaard & Tripic co 1986). The lower 
limit also corresponds to the least massive star on the main se- 
quence for which we predict relatively large surface abundance 
anomalies (see Sect. l5.3l l. All models are shown for the same 
mass loss rate (lO^'^Moyr"'). As will be seen in the next sec- 
tion, this mass loss rate has an effect, though very moderate, on 
the structure via abundance changes (the effect on He abundance 
for instance can increase the depth of the surface convection 

' This question remains unanswered (see iFeldman & Widing|l2003l 
for a complete review). 

Helium is assumed to have a 1 : 1 ratio with high-FIP elements, since 
the observations which suggest th at He would have a ratio o f 1:4 with 
these elements are questioned (see lFeldman & Widing|l2003h . 
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zone). This mass loss rate also allows surface abundance anoma- 
lies which are compatible with observed abundance anomalies 
for Am stars (see Sect. l5.3b . 

The main-sequence lifetime ranges from about 500 Myr for 
the 2.5 Mq model to more than 3 Gyr for the 1.3 Mq model. The 
Fe surface abundance is intimately coupled with the movement 
of the surface convection zone; for the models with the thinnest 
SCZ, Fe is predicted to be overabundant over most of the MS 
lifetime (in the 2.5 Mq model by a factor of about 3). 

All models of at least 1 .5 are marked by a rapid Fe abun- 
dance peak which occurs at the beginning of the main sequence. 
The rise of X(Fe), which is related to the depth change of the sur- 
face convection zone as the star stabilizes (see Fig.|4f), is so rapid 
that it is not resolved in Fig.|4]F. Similarly, there is a sharp spike 
towards the end of the MS for most models which is once again 
correlated with the sharp variation in SCZ depth (in this case the 
most important variation is for the 1.7 model, which has a 
difference of a factor of 2.5 in Fe abundance within less than 
100 Myr). As we will see in the following sections, these varia- 
tions are larger when the mass loss rate is smaller. The HR di- 
agram for these models is shown in Fig.[T] For clarity, only two 
evolutionary tracks are shown from the initially homogeneous 
pre-main-sequence up to the subgiant branch (the others are 
shown from the beginning of the main-sequence, although their 
evolution is calculated from the pre-main-sequence). Diffusion 
and its effects on surface abundances appear well before the ar- 
rival on the main-sequence (a thorough investigation of diffu- 
sion on the pre-main-sequence will be discussed in Vick et al., 
in preparation). 

5.1. Radiative accelerations, internal abundance variations 
and structure 

In Figs. |5] and |6] the radiative accelerations as well as the cor- 
responding spatial abundance variations for a few selected el- 
ements are shown for models of 1.5 Mq and 2.5 Mq. These 
masses approximately coiTespond to the lower and higher Teff 
limits of AmFm stars. 

The MS lifetime of the 1 .5 Mq model is about 2 Gyr whereas 
the 2.5 Mq star has a MS Hfetime of about 520 Myr (see Fig.©. 
The 1.5 Mq model with a mass loss rate of 1 x lO"''*M0yr"' 
(see bottom panel of Fig.jSj was stopped at 575 Myr because the 
solution became numerically unstable. 

Throughout most of the stellar envelope, the mass loss rate 
does not significantly affect the resulting gmd^, as seen in the 
upper row of Fig.|5] However, if mass loss is small enough to 
permit iron peak accumulation below the SCZ, as is the case for 
the 1.50W1E- 14 model, then all gr-^^s will be affected by compe- 
tition in the region where metals have accumulated. 

We present in Sect. 5.1.1 and Sect. 5.1.2 approximate equa- 
tions which will facilitate the interpretation of the detailed solu- 
tions shown in Figs.|5]and|6] which will be discussed in Sect. 5.2. 

5.1 .1 . The interior wind solution 

Consider the approximate solution to Eq. Q, in a regime for 
which the ^ term is small compared to the others (which is true 



over most of the MS lifetime). Then one may write: 
V-[-pDVlnc-Hp(U-HU,,)c] =0 



(18) 



if nuclear terms are negligeable, which is true for the stellar en- 
velope, and 5 w = 0, which is true outside the surface convection 
zone. In one dimension, and in the absence of turbulence (so 



that D becomes D12, the diffusion coefficient in a hydrogen and 
helium background) Eq. ifTSll becomes: 



dc 

-pDi2^ +p(U + UJc 
or 



est 



(19) 



where f^ir) is constant as a function of r and U is approximately 
given by Eq. [9] from .Michaud et al...l983. with a sHght correc- 
tioifl: 
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kT 



1 \ mpg 
2) kT ^ 
d\nT' 



+ aj- 



dr 



(20) 



Here A and Z are the atomic mass and charge number respec- 
tively, and aj is the thermal diffusion coefficient Since the g 
and giad terms are multiplied by A while the abundance gradient 
term (in Eq. |[T9l ) is not, only very large abundance gradients can 
affect the diffusion velocity. Therefore, the internal abundance 
distribution approximately satisfies the relation: 



T{r) = est ^ r^piU + U^)e 



(21) 



which expresses the conservation of particle flux throughout the 
stellar envelope. In this discussion, for the purpose of illustra- 
tion, one may also neglect the term involving because it is 
not dominant in the relevant regions of the envelope (although it 
is always included in our calculations). 

The implications of flux conservation are illustrated in Fig.|7] 
for three atomic species (Ca, Mn and Ni) from the 2.5 Mq model 
of Fig.|6l For all three species, the flux is always positive as may 
be seen from the top row. While giad(Ca) is smaller than gravity 
at some mass shells (middle row), the downward diffusion ve- 
locity of Ca is never larger than the wind velocity (see Fig.fTOb 
so that it is dragged by the wind toward the surface. In general, 
as long as the absolute value of the diffusion velocity is smaller 
than the wind velocity, the concentration simplyincreases in or- 
der to conserve the flux (as implied by Eq. 1211 o This equation 
may then be used here for all three species. As may be seen from 
the bottom row, after 3 Myr, the concentration has adjusted only 
down to AM/M, ^3x10"^. Therefore, below this depth, the local 
flux (top row) mainly reflects the local grad- After about 10 Myr 
however, the abundance has adjusted to carry the flux which ar- 
rives from deep inside the star (down to log AM/M»^ -6). This 
may be evaluated using AM ^ M x t for a mass loss rate of 
lO^'^Moyr"'. After 100 Myr, the concentration has adjusted to 
the flux down to log AM/M, ^ -5 and the concentration becomes 
the mirror image of gi-ad from the surface to that depth as may be 
seen by comparing the middle and bottom rows. There are Ca 
overabundances at log AM/M» - -8.5 and -4.5, which occur at 
giad(Ca) minima (see also Fig.|6]l, while the minimum abundance 
occurs at the maximum of gi-ad(Ca). However, if the wind veloc- 
ity were not larger than the settling velocity, a gradient could de- 
velop to conserve the flux. If this gradient cannot become large 
enough to conserve flux, then the de/dt term (in Eq. IT)), which 
is required to be zero in the kinematic approximation becomes 
important, and the approximation leading to Eq. 11211 becomes 
invalid. 



' ^ had been erronously forgotten in the 2nd and 3rd terms. 

'° This equa t ion is n ot used for the calculations. For more details see 
iTurcotte et al ]| |l998bh . 

" This resembles th e results for oxygen which are discussed in 
iLandstreet elaP ( Il998h . 
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Fig. 4. Evolution of Tetf (a), \ogg (b), luminosity (c), temperature (d) and mass at the bottom of the surface convection zone (e), as 
well as the surface Fe abundance (f) for selected masses that span the range of AmFm stars. All models are shown for a mass loss 
rateof 10 '^Moyr '. The curves are identified in panel (a). Vertically, they have the same order in panels (a), (c) and (f), but are in 
reverse order in panels (b), (d) and (e). 



The net elemental flux toward the surface quickly becomes 
nearly constant in time once the local abundance has adjusted 
itself to conserve the flux. The remaining slow variation in time 
of the flux comes from the variation of grad where the matter 
originated, deep in the star, at AM ^ M x t. 

5.1.2. The kinematic approximation 

By following the analysis of iMichaud & Charlandl (Il986h . the 
surface abundance of elements which are pushed upwards by 
grid and/or dragged by the wind throughout the stellar envelope 
can be approximated by using a simple kinematic equation so 
long as the evolutionary effects and the contribution of |j are 
small ( i.e. there is no significant accumulation). The problem of 
determining surface abundances thereby reduces to a kinematic 
problem. The only required quantity is the local velocity and, 
consequently, the radiative acceleration. Elements which origi- 
nated at ri at f = arrive at r^z, the radius at the bottom of the 
surface convection zone at a time f i given by: 



Because of flux conservation, the flux of an element A entering 
the convection zone at ti is then given by: 

4>{ti,A) = Co(A,ri)(VD,ri + Vw,r,)Pfi-5-, (^3) 

r 

cz 



where co(A,ri) is the initial abundance of A at ri and vdi, Vwi 
and pi are evaluated at ri. This is valid for the region above 
which Vw + vd > and for t < toA, the time required for the 
element A to migrate from the point where + - to the 
bottom of the convection zone. Since diff'usion timescales are 
much longer than convective mixing time scales, the convection 
zone is assumed thoroughly mixed up to the bottom of the wind 
forming region. Within the context of unseparated mass loss, the 
evolution of the abundance of A in the surface convection zone, 
Ccz(A) can then be approximated by: 

Me, = -Mccz(A) + 47Tri4>(ti , A), (24) 

ot 

where Mc, and re- are respectively the mass of the superficial 
convective zone and the radius at the bottom of this convection 
zone. Eq. Il24ll is not used for our calculations, but is used to 
interpret the results. 

5.2. Discussion: Internal variations 

By comparing Figs. |5] and |6] it is clear that the same mass loss 
rate leads to quite different internal concentration variations in 
stars of different mass. There are two main reasons for this. First, 
in the more massive star, the surface convection zone is much 
thinner, therefore the radiative zone is extended upward into re- 
gions where the diffusion timescales are much shorter Second, 
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Fig. 5. [Top row] Radiative accelerations (solid line) and local gravity (dotted line) for a few selected elements in a 1.5 Mq model 
with two different mass loss rates. For the 1.50W1E-13 model, the radiative accelerations are shown for 2 ages (30 and 503 Myr), 
while they are only shown at 500 Myr for the 1.50W1E- 14 model. The vertical lines show the position of the bottom of the surface H- 
He convection zone. The corresponding internal abundance variations at different ages (in Myr) are shown for both the 1.50W1E-13 
model [middle row] and the 1.50W1E-14 model [bottom row]. 



because of the Tejf dependence of the photon flux, most g^^^s 
are stronger in more massive stars. 

For both models, diffusion mainly afi'ects the outer 10 of 
the star's mass; the point above which the eff'ects of atomic diffii- 
sion become visible on Figs.|5]and|6]will be defined as the point 
of separation. Atomic diffusion can also act deeper in the star 
(e.g. in the core), though in A and F main-sequence stars, its ef- 
fects are much smaller below than above the point of separation. 

5.2.1. The 1.5 Mq models 

For the 1.5 Mq model with a mass loss rate of 1 x lO^'-'Moyr"' 
(henceforth designated 1.50W1E-13 for short), internal anoma- 
lies are small. Iron reaches an overabundance of 1.5, and nickel 
of 2.0. Lithium, oxygen and calcium decrease from the point 
of separation up to the bottom of the surface convection zone. 
Dififerences in ^lad lead to diff'erences in their behavior. Due to a 
strong gi-ad immediately below the surface CZ, Ca is most under- 
abundant in that region, since Ca is pushed into the convection 
zone by both g,ad and the wind while little arrives from below 
because of the minimum in its g,-ad at log AM/M, ~ -4.5. At 
500 Myr, its underabundance reaches a factor of 4 while the Li 
underabundance below the SCZ is about -0.3 dex. 

In the lower row of Fig.|5] one sees that when the wind is 
10 times smaller, the internal variations are much stronger since 



the advection by Vwind is not strong enough to prevent elemen- 
tal accumulation. Lithium has an interesting behavior: it has a 
local maximum where its gmd equals gravity. The underabun- 
dances of Li and O drop to -0.95 and -1.35 dex respectively, 
while underabundances greater than 2 dex are reached for S and 
Ca. Through Eq. ET\ . when U^, is smaller, the change in U 
due to giad has a larger effect on concentration. This is partic- 
ularly true for Fe: due to a dip in grad (see Fig.|5]l, it accumulates 
near log AM/M. = -6.5 (where T = 200 000 K). The overabun- 
dances of Fe and Ni reach 1 .25 and 1 .45 dex respectively. The 
implications of this accumulation are analyzed in Fig.|8] The lo- 
cal opacity bump in the region between log AM/M, ^ 6.0 - 6.8 
leads to the appearance of an iron peak convection zone (V^- Vq^ 
changes sign) at around 70 Myr which survives until the end 
of the simulation. Fe and Ni reinforce each other. Only the 
lO^'^^Moyr"' model becomes convective; doubling the mass 
loss rate reduces the Fe accumulation enough for the iron peak 
convection zone not to appear. 

In fact, a noteworthy transition occurs in internal solu- 
tion types between the 1.50W1E-14 and 1.50W2E-14 models 
(Fig.[8]l. In the 1.50W1E-14 model, Fe is overabundant where 
^rad(Fe) > g, whereas Fe is minimal at the grad(Fe) maximum in 
the 1.50W2E-14 model (similarly, Ca is minimal at the giad(Ca) 
maximum in the 2.50W1E-13 model shown in Fig.|7]i. This is 
due to the relationship between v^ind and the settling velocities 
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Fig. 6. [Top row] Radiative accelerations (solid line) and local gravity (dotted line) in a 2.5 Mq model with a mass loss rate of 
lxlO"''^Moyr"'. The corresponding internal abundance variations are shown at three ages [bottom panel]. All curves are identified 
by their age (in Myr). The vertical lines show the position of the bottom of the surface H convection zone. 



for each individual species; if the setthng velocity dominates lo- 
cally, the element will accumulate locally, and the solution will 
behave as does Fe in the 1.50W1E-14 model. When Vwind domi- 
nates throughout the envelope, the solution is determined by flux 
conservation, as shown in Fig.|7] This transition separates the so- 
lution in which the surface abundances reflect matter which is 
advected from deep inside the star (vwind > t^'sett) from the solu- 
tion in which the surface reflects variat ions at the bottom of the 
SCZ (as in the models of I WatsonllT97Tl and I Alecianll 1 996h . This 
transition occurs only for elements whose ^rad has a minimum 
smaller than gravity close to the bottom of the SCZ. 

A large enough mass loss rate also affects the position of 
the bottom of the surface H-He convection zone by keeping He 
from sinking. Since it also modifies the accumulation of metals, 
mainly iron and nickel, its effect is complex: the relative position 
of the bottom of the convection zone for the three mass loss rates 
in Fig.lHis different at 70 and 500 Myr. 

The inversion of the molecular weight gradient (bottom row 
of Fig. [8] at 500 Myr), which eventually follows the appear- 
ance of the iron accumulation around 200 000 K, has been sug- 
gested to affect the presence of the iron peak convection zone 
dTheado et al.l l2009h . In the present calculations however, the 
iron peak convection zone appears (as seen at 70 Myr) before the 
inversion appears so that while the size of the convection zone 
could be affected by the jj gradient inversion, its appearance can- 
not be affected. This will be further discussed in Sect.|8] 

5.2.2. The 2.5 Mq models 

Large overabundances of iron peak elements are obtained in the 
2.5 Mq models (Fig.|6]i: Ni and Fe are overabundant by factors of 
7.5 and 4 respectively. Calcium on the other hand is found to be 
up to 10 times underabundant around log AM/M* =-6.5, where 
its grad is near a maximum. This may seem counter intuitive but 
is a consequence of Eq. 1211 . Since the flux is conserved, where 
vd is large and positive, c decreases, which is precisely what 



happens for Ca just below the convection zone (see Figs. |6] and 
|2]and the last two paragraphs of Sect. l5.1.T] ). 

Wind and diffusion velocities are compared in Fig. [10] The 
wind velocity decreases more rapidly than the diffusion veloc- 
ity as log AM/M, increases. For instance, the absolute value of 
the "^He settling velocity is 50 times smaller than the larger of 
the two wind velocities at logAM/M» = -12 but equals it at 
log AM/Mt = -4.5. For smaller mass loss rates the settling ve- 
locity dominates closer to the surface. Calcium has an upward 
diffusion velocity over the interval -9 < log AM/M, < -6 and 
it is up to ten times larger than the wind velocity. Fe and Li have 
smaller upward diffusion. Those velocities determine local con- 
centration via Eq. 1211 (see the end of Sect. l5.1.TI ). 

Fig.|9] allows an analysis of the effect of varying mass loss 
rates on surface convection zones and to distinguish the effects 
of mass loss from those of turbulence (see Sect.|6]l. The abun- 
dances and giad of both Fe and Ni are presented since they are 
the main contributors to the appearance of the iron peak convec- 
tion zones in 1.5 Mq models. In 2.5 Mq models however, iron 
convection zones do not appear for the two mass loss rates con- 
sidered. The Fe and Ni abundances remain slightly below solar 
where grad(Fe) and grad(Ni) are largest so that they do not lead 
to an increase of opacity where they could contribute most to 
opacity. The wind velocity is ten times larger than the diffusion 
velocity of Fe for log AM/M, > -8 so that Eq. 1211 forces the 
solution to be nearly constant over that range. This is to be con- 
trasted with their behavior in 1.5 Mq models. 

Even a small mass loss rate difference can have important 
effects on the stellar structure in 2.5 Mq models. For both mass 
loss rates considered, the He I convection zone disappears early 
in the evolution, at around 30 Myr (see the fifth row of Fig.|9]l. 
The Hell convection zone on the other hand only disappears 
in the model with the smaller mass loss rate, the 2.50W5E-14 
model, at around 200 Myr. This is because the inward diffusion 
velocity of He dominates the wind velocity at a shallower depth 
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Fig. 7. Comparison of the normalized local flux with radiative accelerations and internal abundances for 3 elements at 4 different 
ages (in Myr) for a 2.5 Mq model with a mass loss rate of 10"'^ Mgyr"' . 



than in the 2.50W1E-13 model (see Fig.fTOli. so that a more pro- 
nounced He underabundance develops to conserve the flux. 



5.3. Surface abundance variations 
5.3.1 . The 1 .5 Mq models 

The evolution of surface abundances is shown for several 1 .5 Mq 
models in the upper two rows of Fig.[TT] For the 1.50W1E-12 
model, the largest anomaly encountered, which is for Ni, reaches 
a mere 0.05 dex, and it is quickly flattened after a brief period. 
Such a mass loss rate effectively wipes out any surface effects of 
chemical separation. In the 1.5 W IE- 14 model, some of the cal- 
culated underabundances are very large (greater than 2 dex for S 
as well as C and Si, which are not shown). With the exception of 
Fe, all iron peak elements are overabundant throughout the sim- 
ulation. Both Ni and Mn become more than 1 dex overabundant 
at around 500 Myr and 50 Myr respectively. 

Iron surface abundances are particularly interesting for this 
model since it is the only one to have an iron surface underabun- 



dance which spans 

500My]B 

This is due to the accumulation of 
iron under the surface convection zone, and the subsequent ap- 
pearance of an iron peak convection zone (see Fig. [8]). Although 
not shown here, a separate iron peak convection zone is expected 
for all simulations for masses between 1.47 and 3.0 Mq 
when the mass loss rate is equal to or below ~ lO^'^^Moyr"'. 

The 1.50W2E-14 and 1.50W5E-14 models generate sur- 
face iron overabundances of about 0.5 and 0.35 respectively at 
300 Myr. The calcium abundance is particularly interesting since 
for all models, an overabundance is predicted at the beginning 
of the main sequence evolution followed by an underabundance. 
The underabundance is present over a much larger fraction of 
the evolution than the overabundance. The Ca calculations ar e 
compared in Fig. [12] with those carried out by lAlecianI (Il996l) . 
His calculations were undertaken in static stellar models; there- 
fore, evolutionary effects were not included. There are also slight 
difference in Tefj between our calculations (Teff ^ 7200 K at 
the beginning of diffiision but decreases as the star evolves) and 



We may assume however that if the simulation had not encoun- 
tered significant instabilities, the shown iron abundance trend would 
have continued and iron would have become overabundant. 
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Fig. 8. Internal abundances (Fe, Ni), radiative accelerations (Fe, Ni), Rosseland opacity, the difference between the radiative and 
adiabatic temperature gradients as well as the mean molecular weight per nucleus for three 1 .5 Mq models with different mass loss 
rates (10"''*MQyr"' [dashed line], 2 x lO^'^^Moyr"' [dot-dashed Hne] and lO^'^Moyr"' [solid hne]) at 70Myr (left column) and 
500 Myr (right column). At this mass, the Fe peak convection zone (between log AM/M, = -6 and -7) only appears for the mass 
loss rate of 10 '"^Moyr The H-He convection zones remain linked for all mass loss rates. 



his simulations (reff - 7500 K). There are also differences in 
giad(Ca): his calculated ^iad(Ca) has a peak which is 10 times 
smaller than obtained in a model of similar T^ff with our code. 
Results are compared in Fig. [12] for three mass loss rates. In all 
three cases, the maximum anomaly as well as the overall be- 
havior of the curves correspond well. The overabundance peaks 
are however not quite so wide in our calculations as in his. The 
agreement seems satisfactory. 



5.3.2. The 2.5 Mq models 

In the 2.5 Mq models (lower two rows of Fig.fTTI). the surface 
abundances of heavier metals (say Z > 17) can be character- 
ized by two distinct episodes: a steep abundance spike, followed 
by a smooth decline. The initial peaked episode around 10 Myr 



shows matter which was above log AM/M» ^ -6 at the onset of 
diffusion, advected to the surface by the stellar wind. 

As discussed at the end of Sect. l5.1.Tl the initial flux dis- 
tribution reflects giad since, originally, the concentration is the 
same throughout the envelope; then, as time elapses, the internal 
abundances naturally evolve in such a way that flux is conserved 
throughout the envelope (see Fig.|7]i. In so far as the dc/dt term is 
negligible (i.e. flux is conserved), the surface reflects the point in 
the initial flux distribution which is dragged by v„ind. Therefore, 
when the initial internal flux profile of an element has important 
variations near the surface due to strong variations in (as is 
the case for Ca and Ni, see Fig.|6]and|7]), these variations appear 
on the surface in a time which is related to the mass loss rate. 

Quantatively, in the 2.50W1E-13 model (solid line), a given 
element's surface abundance at 100 Myr depends on the initial 



14 



M. Vick, G. Michaud, J. Richer, and O. Richard: AmFm and lithium gap stars 



41 



> 



.1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I4J 



Fe 



I I I I 



t 1 " " 


" " 1 1 


Fe J 




\ 






'■\^ \^ 




T 1 . I I I 


1 1 1 1 1 1 " 


1 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I l-U 



Ni 



.-,\ 



\ 



t I I I I 



1 

'- -\ - 


' " 1 1 


1 1 1 1 1 1 1 1 


' 1 


' 1 ' 

Ni J 












1 


. . . 1 1 






. 1 




1 


1 1 




1 ' 








Ni : 


'Jr' 




log « 




\l < 


1 1 


1 


1 ; 



i 


1 


' " 1 J-'J ' ' 


' 

. Ni : 










i , i, 


1 


1 1 1 1 1 


T 



1 

0.1 
0.0 
-0.1 
-0.2 





1 1 1 1 1 1 


1 


1 1 " 


1 1 1 1 










































Till 





1 


1 >^ 





1 " 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 M 

A 


r, 


\ ^ 



1.3 
f 1.2 
1.1 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



I I I I I I I 



1 1 ' " 


1 1 1 ' " 


1 1 1 




















T 1 . . . 


1 1 1 r=tai 


1 1 1 




r "\\ 1 ' " 


1 1 1 ' " 


1 l-l 


^\ 






- ^..^ 

- 1 1 






-H He I 

1 1 . . . 


He II Fe'.Ni^^X./^''""'^ 

1 1 1 1 1 1 


1 1 1 




V 1 ' " 


1 1 1 ' " 

.■■/ 


1 1 1 















r 1 , , , 


1 1 1 1 1 1 


1 1 1 


12 -10 


-8 -6 -4 





-12 -10 -8 -6 -4 

log AM/M. 

Fig. 9. Internal abundances (Fe, Ni), radiative accelerations (Fe, Ni), Rosseland opacity, the difference between the radiative and 
adiabatic temperature gradients as well as the mean molecular weight per nucleus for two 2.5 models with different mass loss 
rates (lO^'^Moyr"' [dotted line] and 5 x lO^'^^Moyr"' [dashed line]) as well as a model with turbulence [solid line] at 30Myr (left 
column) and around 250 Myr (right column). The maxima in the Vf-Vad row are mainly due to opacities from H, He I, Hell and 
Fe/Ni from left to right respectively. 



7 -2.0 P 
w 

6-4.0 

o 

^ -6.0 



O -8.0 



I I I I I I II I II I II I II 



"■\\ 



"l I I I I I I ' 



-9 -6 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 

^Li - 














1 1 1 1 


III 




-9 


-6 



1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 




- 




V\ 




'■\ \ - 










"i 1 1 


1 1 1 



II 1 II 1 


1 1 1 1 1 1 1 1 1 1 1 


::^^\^ 


Ca - 




\\\ ~ 




\ A 






"i 1 1 1 1 1 


1 1' " 



J 1 1 , 1 


1 1 1 1 1 1 1 1 1 1 1 




Fe - 














"i 1 1 1 1 






-9 -6 -9 -6 

log AM/M, 



-9 -6 



-9 -6 



Fig, 10. Wind velocities (long dashed line: lO^'^Moyr"'; dashed line: 10^'^Moyr"') and diffusion velocities (solid when positive, 
towards the surface, and dotted when negative) of a few selected elements in two 2.50 models at 300 Myr. For most species, 
wind velocities decrease more rapidly inwards than diffusion velocities. 



M. Vick, G. Michaud, J. Richer, and O. Richard: AmFm and lithium gap stars 15 




-3 -2 -1 -3 -2 -1 -3 -2 -1 -3 -2 -1 



log Age(Gyr), M.=2.5 

Fig. 11. Surface abundance variations for five 1.50 Mq models [top 2 rows] as well as three 2.50 Mq [lower 2 rows] models with 
Zq - 0.02 and varying amounts of mass loss (in Moyr"^.). The model was ended at around 550 Myr due to numerical instabilities 
which were caused by the severe underabundances of many elements (of which S is shown). It is also the only model which leads 
to the appearance of a separate iron peak convection zone. In the bottom panel, the two models with smaller mass loss rates lead to 
a surface metalUcity which has a maximum around 10 Myr before falling below Zq around 100 Myr. 



flux variations (caused by g^^d variations at f = OQ at a depth 
of AM ^ M X f = lO^'^Moyr-i x 10** yr ^ IO^^Mq. In Fig.[n] 
the slight bump in surface Ni abundance around 100 Myr for the 
2.50W1E-13 model reflects the small bump in the initial flux 
distribution of Ni around log AM/M, ^ -5 (see Fig.|7]l. 

The previous example is an application of Eq. Il22l and shows 
tiiat the time it takes for a given internal variation to reach the 
surface is inversely proportional to the mass loss rate. This is 
further illustrated by the fact that the surface Ni abundance peak 
around 10 Myr appears earlier as the mass loss rate increases (see 
Fig.fTTTi. The dilution by the convection zone (Eq. Il24l ) is rela- 
tively small because it has a relatively small mass. Also, the rapid 
variations seen around this Ni abundance peak reflect the many 
internal flux variations which were above logAM/M, ^ -6 at 
f = 0. For a mass loss rate of lO^'-^Moyr"', these variations all 
reach the surface within 10^ yr. 



Remember that at / = 0, the flux distribution is essentially pro- 
portional to the local abundance multiplied by the local g^^^ since the 
composition is homogeneous. 



As a consequence, the nearer an internal flux variation is to 
the surface at f = 0, the quicker it appears at the surface, and the 
quicker it disappears. A flux variation that spatially spans from 
the surface to log AM/M» = -7 at f = appears and disappears 
at the surface in less than one million years for a mass loss rate 
of lO^'^'Moyr-i. 

The amplitude of the surface variations also depends on the 
mass loss rate; as the mass loss rate increases, each internal vari- 
ation of the flux at f = manifests itself at the stellar surface 
with a smaller amplitude since more enriched/depleted matter 
is evacuated from the SCZ (Eq. Il24ll ). This is illustrated by the 
decreasing amplitude of the initial surface abundance spikes for 
Mn and Ni as the mass loss rate increases (Fig.fTTI). 

The conservation of the flux down to log AM/M» ^ -5.5 to 
-6.5 (depending on the age and mass loss rate) has another im- 
portant consequence: diffusion between the surface convection 
zones (i.e. between the H and either of the two He convection 
zones) has practically no influence on the surface abundances. 
Matter which originates from between the H and He convec- 
tion zones (log AM/M, ^ -10) appears at the surface within lO-' 
years for a mass loss rate of lO^'^Moyr"'. The detailed calcu- 
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Fig. 12. Evolution of Ca surface abundances for 1.5 Mp mod- 
els wit h different mass loss rates. The calculations from lAlecianI 
119961) are shown in bold. 
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Fig. 14. Surface abundance anomalies at five ages (5, 9, 70, 
250 and 524 Myr) for a 2.50 Mq star with a mass loss rate of 
10^'^Moyr-'. 



lation of chemical transport between these zones is then not re- 
quired in order to accurately obtain the surface abundance solu- 
tion. To verify the accuracy of this assertion a model was calcu- 
lated with homogenized abundances between all surface convec- 
tion zones, and it was found that the surface solution was practi- 
cally identical to the surface solution obtained when separation 
was allowed between SCZs. After a mere 3 Myr of a star's main- 
sequence life with a mass loss rate of lO"'-' Moyr"', anomalies 
which appear at the surface reflect the separation which occurs 
below logAM/Mt — -7. Hence, for most of the main-sequence 
lifetime of models compatible with observations, the surface 
abundance solution depends on the separation which takes place 
around log AM/Mt ^ -5.5 to -6.5. 

Th is is analogous to the turbulence models of iRicher et alJ 
( I2OOOI) . in which surface abundances depend solely on the sepa- 
ration which occurs below 200 000 K. In our calculations, how- 
ever, abundance variations are present throughout the stellar en- 
velope because no mixing is enforced outside of convection 
zones. In that respect, in the presence of mass loss, it is clear that 
the chemical separation responsible for the AmFm phenomenon 
involves up to log AM/M* = -5 of the star's mass. 

5.3.3. The effect of Z, age and Teff 

In Fig. [131 one sees that the main features of the time evolution 
of surface abundances are similar for three different initial metal- 
licities. However varying the initial metallicity can have an im- 
portant effect on the amplitude of surface anomalies, though not 
for all elements. For elements such as CNO, varying Zq has rela- 
tively little effect on the surface abundance anomalies since these 
elements are not supported by the radiation field. For heavier 
iron peak elements, supported by the radiation field, the ramifi- 
cations are much more apparent since their lines are often sat- 
urated and flux sharing becomes prevalent. In fact, the lower 
metallicity models have larger anomalies, when compared to the 
original abundances of the model, since the same flux is shared 
among fewer atoms. Also, the abundance peaks appear earlier 



during evolution as Z decreases since the g^^d profiles with re- 
spect to AM/Mt are shifted toward the surface as Z decreases. 
However the situation is different if one compares to a fixed set 
of abundances. By comparing Figs.[TT]and[T3]one notes that a 
reduction of mass loss by a factor of 1.5 approximately com- 
pensates for a reduction of Zq from 0.03 to 0.02 for the absolute 
final abundance of Fe and other elements which are supported. It 
however amplifies the effect of the reduction for elements which 
are not supported by ^lad such as CNO. 

Figure[T4] shows the surface abundances as a function of 
atomic number for a 2.5 Mq model at five different ages. A 
large overabundance of Ti ( 1 . 3 dex) occurs as early as 5 Myr be- 
cause grad(Ti) has a maximum near the surface, and the bump it 
causes in the original flux is dragged by the wind. The same 
occurs for Cr and Mn but to a lesser extent. There is also a 
slight overabundance of Li which appears at this age. At 9 Myr, 
the Ti overabundance has weakened and most iron peak ele- 
ments reach their highest values: Fe peaks at 0.65 dex and Ni 
peaks at 1.0 dex. Throughout the star's main-sequence evolu- 
tion, the surface abundances change only slightly, as highlighted 
by the similarity in the curves at 70 and 250 Myr. Until the 
final dredge-up (i.e. the post turn-off increase in SCZ mass), 
which begins around 520 Myr, the iron peak abundances de- 
crease only slightly. For instance, Fe only drops 0.15 dex over 
the entire main-sequence lifetime. The changes will be slightly 
more significant for underabundant elements because they are 
not supported by the radiation field, and so the amplitude of 
the underabundances is largely dependant on time (due to set- 
tling). For CNO, the surface underabundances appear around 
9 Myr and continuously decrease down to about -0.40, -0.38 
and -0.35 dex for C, N and O respectively, after which dredge- 
up begins. At all times, the surface abundances reflect the matter 
which is exposed by the wind at the specific age. With that in 
mind, one notes from the curve at 524 Myr, that the material at 
log AM/M, ^ -2 has not yet reached the surface since LiBeB 
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Fig. 13. The effect of varying initial metallicity on the evolution of abundance anomalies at the surface of a 2.0 Mq star with a mass 
loss rate of lO^'-'Moyr"'. Models are shown for Zq - 0.01 (solid), Zq - 0.02 (dashed line) and Zq - 0.03 (long dashed line). The 
original abundance of Li is assumed independent of Zq. 



abundances have not yet dipped due to material advecting to the 
surface from regions where they burn (see Fig.|6]l. At this age, the 
surface convection zone extends down to log AM/M* = -4.5, 
and while iron peak elements are now underabundant, an over- 
abundance has appeared for S (and a few species of similar 
atomic mass) because of the dredge up of the maximum of its 

abundance at logAM/M« 5.2 (see Fig.|6] one also notes the 

underabundances of Fe and Ni near log AM/M» ~ -5 which 
explain their underabundance at 524 Myr). This time dependent 
behavior is compared to models with turbulence in Sect.|6] 



Figure[T5] shows the superficial abundances at 500 Myr 
for models of different masses with the same mass loss rate 
(5 X lO"''*M0yr"' ). The 1.3 Mq has much smaller anomalies 
than the others. At a given age, stars of 1.5 to 2.3 Mq have very 
similar underabundances of elements from He to Ne, since for a 
given mass loss rate, they essentially depend on time. Iron peak 
overabundances increase as Teff increases. This is partly caused 
by the increase of g^ai with T^ff as well as by the reduction of the 
mass of the surface convection (if the SCZ is more massive, the 
anomalies will be reduced by dilution). This behavior is differ- 
ent from what is obtained in turbulence models, and will also be 
discussed in Sect.|6] 
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Fig. 15. The effect of varying Tetf (or stellar mass) on the surface 
abundance profiles of 5 models at 500 Myr for a mass loss rate 
of 5 X lO-i^Moyr-i. 



5.3.4. Stars of the lithium gap 

The surface Li abundances for stars of the lithium gap without 
mass loss are shown in Fig. [16] The models heavier than 1 .45 Mq 
ceased converging before the end of their main-sequence life- 
time because of numerical instabilities. In the absence of mass 
loss, the H-He and iron peak convection zones are split by the 
appearance of a thin radiative zone for all models of 1.47 Mq 
or more. Variations within this layer appear at the surface of the 
1 .47 Mq model near 200 Myr. 

A close inspection of Fig.[T6lshows a disctinct separation be- 
tween the 1.46 Mq and 1.47 Mq models; the 1.46 Mq model be- 
haves more like the 1 .45 Mq model while the 1 .47 Mq model's 
behaviour most resembles the 1.48 Mq model. Again, this is be- 
cause of the appearance of a radiative zone immediately above 
the iron peak CZ in the 2 heavier models. 

The abundances obtained near the age of the Hyades 
open cluster ranged from -0.25 dex for the 1.37 Mq model to 
-1.25 dex for the 1.46 Mq model. Lithium is underabundant for 
all models throughout evolution. 

By comparing these results with those shown in Fig. 6 of 
[Richer & MichaudI (Il993h . one first notices that the curves have 
very similar behavior in time; the minima occur at nearly the 
same age and the curve shapes are nearly identical. However, 
the underabundances obtained in our calculations are systemati- 
cally smaller. For instance, their 1 .43 Mq model is 250 times un- 
derabundant around 1 .23 Gyr, while our 1 .43 Mq model is only 
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14 times underabundant at the same age. A careful analysis of 
the results determined that the difference was principally due to 
the difference in the mass of the SCZ. If the convection zone is 
homogeneous and emptied through its bottom, then: 



10.000 



dt 



Xwip 



(25) 



where X is the mass fraction in the SCZ (and on the surface), 
AMcz the mass in the SCZ, w, the diffusion velocity and p the 
density immediately below the SCZ. By integrating this equa- 
tion and comparing diffusion time scales for the different values 
of Mcz, the differences in surface lithium abundances for models 
of the same mass can be explained by the difference in M^z- This 
difference in CZ mass is mainly due to the inclusion of metal dif- 
fusion in our models. In fact, the increased opacity in our models 
due to heavy metal accumulation considerably increases opacity 
at the BSCZ, thereby extending it inwards. 

In Fig. [17] the surface lithium abundances are shown for dif- 
ferent isochrones of the models shown in Fig. [16] The gap posi- 
tion evolves toward cooler temperatures with time. At lOOMyr, 
only stars of 6900 K have significant underabundances, while 
at 625 Myr, stars as cool as 6700 K have important underabun- 
dances. 

At 100, 200 and 400 Myr, the gap obtained in the present 
calculations is 50-80 K hotter than the gap of the same metal- 
licity shown in Fig. 8 of [Richer & MichaudI (Il993h . Our gap is 
also slightly deeper, as the surface lithium on the hotter end of 
our isochrone s keeps d ecreasing where surface lithium in the 
isochrones of lRicher & M ichaudI d 19931) start increasing. This 
is a result of the smaller grad(Li) obtained in our models (see 
Sect. l3.1.2] in particular Fig.[3]i: the competition with Fe de- 
creases ^i-ad(Li) by a factor of 2, so that it doesn't reach g in 
the hotter models of our isochrones. Therefore, lithium is not 
supported in our diffusion only model. 
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Fig. 16. Evolution of ^Li surface abundances for models with 
masses ranging from 1.37 Mq to 1.48 Mq with no mass loss 
and Zo - 0.02. The vertical line indicates 625 Myr, the age 
of the Hya des open cluster. This can be compared to Fig. 6 of 
[Richer & MichaudI (119931) . 
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Fig. 17. The T^ff-Xu relation at different ages for the models 
shown in Fig. [16] Each point represents a calcu lated model. This 
plot can be compared directly to Fig. 8 of iRicher & MichaudI 
(1993) . 



5.3.5. Separated winds 

The effects of various separated wind configurations on surface 
abundances of 1.4 Mq models are illustrated in Fig. [18] Along 
with the 1.40W1E-14 unseparated mass loss model and the pure 
diffusion model (no wind), three cases of separated winds were 
considered, as described in Sect. 14. 31 

In constrast to unseparated mass loss (dotted line) which, 
when comparing to the pure diffusion model (solid line) in 
Fig. 1181 reduces surface anomalies for all elements, separated 
mass loss can affect underabundant and overabundant elements 
differently. In case 1, for which all metals are ejected with the 
same relative concentration as in the atmosphere whereas H and 
He remain bound, a mass loss rate of 1 x lO-"'Mo/y0 leads 
to smaller overabundances and larger underabundances than 
unseparated mass loss. This is because in case 1, the internal 
wind term this mass loss leads to is negligeable (see Eq. [6]l. 
Therefore, for underabundant elements, the surface convection 
zone is drained from its bottom (because the wind is too weak 
to support downward diffusing elements), as well as from its top 
(i.e. mass loss at the stellar surface). Overabundant elements, on 
the other hand, are depleted through the surface, while the bot- 
tom of the surface convection zone is replenished by atomic dif- 
fusion. When the mass loss rate is increased to 5 x lO^'^Moyr"' 
for case 1, the overabundances can rapidly evolve into under- 
abundances since, without the wind providing sufficient replen- 
ishement from deeper inside the star, the elemental depletion 
at the surface quickly dominates the replenishement at the bot- 
tom of the convection zone. For the model with a mass loss 
rate of 5 x lO^'^Mgyr"' underabundances of Li and Fe reach 
-1.5 dex and -0.95 dex respectively at 570 Myr, whereas for the 
lO^'^Moyr"', the Li underabundance reaches -0.55 dex around 
625 Myr, while Fe has an overabundane of 0.45 dex. Generalized 
underabundances are attainable in the context of case 1. 

In case 2, all low-FIP elements are depleted 4 times more 
rapidly than H and other high-FIP elements. Consequently, 
in comparison to the anomalies obtained for the 1.40W1E-14 
model, the anomalies for the low-FIP elements (of which Li, 



This is smaller than the rate of metal depletion for the unseparated 
model 1.40W1E-14, sinceZxM = 1.999 x lO^'^Moyr-'. 
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Fig. 18. Surface abundance evolution for selected elements of 1.40 M© models with different types of wind solutions: one model 
with no mass loss (solid line), one model with a homogenous mass loss rate of l.E-14MQyr"'(dotted line), and 4 models with 
separated mass loss. For the separated winds, the configuration is indicated by a number (1, 2 or 3) which is explained in Sect. l4.3l 
and in the text below. 



Mg, Ca, Mn, Fe and Ni are shown) are greater for underabunt el- 
ements and smaller for overabundant elements, as in case 1 . The 
Fe overabundance is reduced to 0.2 dex around 625 Myr, and the 
Li underabundance reaches -0.6 dex at the same age. However, 
for O, which is a high-FIP element, the curve is tucked in be- 
tween the pure diffusion and unseparated mass loss curves (and 
is therefore impossible to see on the figure), which is what could 
be expected, since the rate of depletion as seen by O, a high-FIP 
element, is smaller than lO^^'^Moyr"' because the lost mass has 
a higher concentration of low-FIP elements (see Sect. 14. 31 ). 

Finally, in case 3, for which H is included among the low- 
FIP elements, the results are quite different: the depletion of H 
at the surface leads to an important increase of He. Since the He 
abundance at the surface remains slightly above its original value 
throughout the simulation, the stellar structure is changed as the 
H-He SCZ is much deeper on the main-sequence than for any of 
the other models. Around 625 Myr, Fe is barely underabundant, 
while Li is underabundant by -0.5 de?|3- At 625 Myr, the model 
in case 3 is 1 00 K hotter than the three other models at the same 
age. 

6. Comparison to turbulence models 

As previously mentioned, models with turbuleri ce (iRicher et all 
l200aiRicharderan.l200ll:rMichaud et al.Ll2005[) have been quite 
successful at reproducing observed properties of AmFm stars 
both on the surface (i.e. abundance anomalies) and in the in- 
terior (pulsation properties of 6 Scuti stars, see iTurcotte et alj 
I2OOOI) . These models with turbulence put forth a scenario in 
which chemical separation takes place below 200 000 K, as op- 
posed to the scenario in which chemical separation immediately 
below the H con vection zone i s resporisible for AmFm abun- 
dance anomalies ( Watsonlfl 97 it I Alecianll 1 9961) . The models cal- 
culated with mass loss offer an alternative scenario which is in 
fact a hybrid of both these scenarios. On the one hand, no exter- 
nal mixing is enforced outside of convection zones, which allows 
for chemical separation to occur throughout the stellar interior 

For a given element, [X/H] is about 0. 1 dex larger than log(Xs„jf/A'o) 
because of the surface underabundance of H. 



However, for most of the main-sequence lifetime, the surface 
abundance is modulated by matter which is advected from deep 
inside the star (while in turbulence models, it is mixed to that 
depth). In order to determine which of these scenarios must be 
favored, it is paramount to differentiate the models so as to en- 
able observational tests. 

Even when surface abundances are quite similar in turbu- 
lent and mass loss models, the interior behaves differently. This 
can be seen by comparing the 2.50W1E-13 and 2.50T5.2D1M-4 
models in Fig. [19] and Fig.|20l In Fig. [19] which compares the 
internal concentrations of all 28 elements for a model with tur- 
bulent mixing and a model with mass loss, there is a stark con- 
trast in the internal distribution of elements heavier than Ar, 
even though the surface abundances are quite similar for most 
atomic species (see Fig.l20b. The Fe abundance in the inter- 
val log AM/M, ^ -6 to -5 is 4 - 5 times larger in the tur- 
bulence model. The same can be said for Ni which, between 
logAM/M» - -7 and -5, differ by about a factor of 10. This 
has an important effect on the local Rosseland opacity. These 
internal differences could therefore allow for some asteroseis- 
mological tests which could help differentiate between the two 
models. 

With respect to pulsations, one note s that around 250 Myr, 
the age of Sirius (iLiebert et al ] I2005h . the Hell convection 
zone is still present in the turbulence model as well as in the 
2.50W1E-13 model, but not in the 2.50W5E-14 model (see 
Fig.|9]l. As seen in the log kr panel, opacities are quite simi- 
lar in the Hell convection zones for the 2.50W1E-13 and the 
turbulence models, therefore, the mass loss model could be 
compatible with ob served pulsation properties of 6 Scuti stars 
(iTurcotte et al.ll2060 ). Further investigation is required to deter- 
mine if the He II opacity bump in the 2.50W5E-14 model is suf- 
ficient to drive kappa mechanism oscillations. 

During evolution there are phases in which important sur- 
face abundance differences differentiate models with mass loss 
from those with turbulence. For instance. Fig. [14] can be com- 
pared directly to Fig. 15 of Richer et al. (2000). The most glar- 
ing difference is the surface behavior shortly after the onset of 
diffusion. In the turbulence model, abundance anomalies appear 
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Fig. 19. Gray coded concentrations of two 2.50 Mq models at 500 Myr. The scale of the radius is linear although the value on the 
left hand side indicates log AM/M,. The contrast is identified in the right insert and spans the interval from 0.3 to 1.7 times the 
original elemental abundances. The anomalies affect the outer 25 % of the radius. Nuclear effects appear for lighter elements near 
the stellar core. For the calculations shown in the left panel, the competing process is mass loss while it is turbulence in the right 
panel. For these two models, the surface abundances are nearly the same (see Fig.l20li. 



slowly and gradually at the surface, whereas large anomalies ap- 
pear at the surface as early as 5 Myr in the mass loss model. This 
is because all chemical separation that occurs near the surface, 
where timescales are short, will rapidly be advected to the sur- 
face by the win43 (see discussion in Sect. l5.3.2] i. In the turbu- 
lence models, mixing is enforced throughout the upper envelope, 
and effectively prevents any chemical separation near the sur- 
face. Likewise, the abundances immediately following the turn- 
off are also different; while overabundant iron peak elements 
such as Mn and Fe can become underabundant as the SCZ ex- 
poses regions in the envelope with important gradients (as is the 
case at 524 Myr), the internal variations in the turbulence context 
are much less important (as an example, compare the Fe under- 
abundanc e around log AM/M, ^ -5 in our Fig. |6] and Fig. 4 of 
[Richer et al. 2000), and so the variations at the surface will also 
be smaller. On the main-sequence however, the differences be- 
tween the two models are relatively quite small. 

Similarly, for the T^m depende nce. Fig. [15] can be compared 
directly to Fig. 16 of ' Richer et alj ( 12000) . The most glaring dif- 
ference is the behaviour of iron peak elements. By comparing 
the models with masses between 1.7 and 2.3 Mq in both figures, 
one notes that in the turbulence regime, the iron peak surface 
distribution is the same for all TeffS at a given age, while there 
is more significant variations for the same elements in the mass 
loss models within the same stellar mass interval. On the con- 
trary, the lighter elements show less variations in the mass loss 
models in comparison to the models with turbulence in the same 
mass interval. 



7. Comparison to observations 

In order to constrain stellar models, and to determine whether 
turbulence or mass loss is the dominant macroscopic process re- 
ducing surface abundance anomalies, it is imperative to compare 
our results with observations. In order to carry out an accurate 
comparison one needs to constrain age, mass and initial compo- 
sition, which as shown in the previous section, all affect surface 



This will be discussed in further detail in a forthcoming paper 
which will look into the effects of diffusion on the pre-main-sequence. 



abundances. To reduce the arbitrariness of the comparison, we 
chose three open cluster stars for which we have a good evalu- 
ation of the initial metal content as well as of the approximate 
age. We will also compare our results to the field star Sirius A 
and the binary system oLeonis. 

In the following sections, [N/H] has its usual meaning: 

[N/H] = log(N/H), - log(N/H)o. (26) 

As mentioned in S ect.|2l this paper is pa rt of a series of pa- 
pers starting with iTurcotte et al.l ( Il998bl) . where the mixing 
length used was calibrated using the Sun for given bound- 
ary conditions, as well as helium and metal content. For con- 
sistency, the same boundary condition, solar composition and 
mixing length are used for Pop I stars. Our models lead to 
abundance variations, or anomalies, relative to those original 
abundances. A number of observers have similarly determined 
anomalies by differential methods with respect to solar abun- 
dances. Moreover, in the model atmospheres used for abundance 
determina tions, most observer s used t he solar abundance mix ei - 
ther from lAnders & Grevessd (11989^ or 'Greve sse etaP (Il996l) . 
Because determinations were sometimes obtained with different 
solar photospheric abundances by different observers and some- 
times with differential methods, their abundance determinations 
relative to the Sun (i.e. the anomalies) are used, when available, 
rather than absolute abundances. The uncertainty that inaccura- 
cies in solar abundances lead to will be discussed in Sect.[8] 

7.1. Field Stars 

Sirius A is the most studied hot Am star (reff - 9800 K), with 
a mass of about 2.14 Mq and an age of approximately 2 50 Myr 
(Gatewood & Gatewood 1978; see also Sect. 4.1 of Riche r et aH 
2000). Figure|20]shows surface abundance determinations for 19 
chemical species (16 of which are included in our calculations) 
from 8 different papers. 

First, for most elements, there is considerable scatter among 
observers. For instance, there is a 0.3 dex difference in Si 
abundance, as well as a 0.4 dex difference in Fe abundance, 
which is the most carefully determined element. We compared 
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Fig. 20. Observations of the surface abundances of SiriusA 
(also known as a CMa, HR 2491 or HD 48915). Circles, 
iSurkhart & Coupry (199^); Upright ope n triangles, 
Robv & Lambert! (|199Q|) : inv erted op en trian gles, \ha.rabex\. et alj 



1982h : thr ee-point stars, iLemkd (1 19891); inv erted three- 



point stars, [L emke (1990); squares. Hill (1995); diamonds, 
iHill & Landstreet ( 1993 ); asterisks, Hui-Bon-Hoa et al. ( 199?) . 
Calculated values are shown for 3 models with mass loss 
(2.50W5E-14, 2.50W1E-13 an 2.501E-12), as well as the 
model with turbulence ( 2.50T5.2D1M-4) w hich represents the 
best fit from Fig. 18 of [Richer et al.l (l2000h . All models were 
calculated with a solar (Z = 0.02) initial metallicity and are 
shown at an age of 250 Myr Vanadium, scandium and cobalt 
are grayed out since they are not included in our calculations. 
The internal abundance distributions of the 2.50W1E-13 and 
2.50T5.2D1M-4 models are shown in Fig.[T9l 



this data to 3 models with mass loss (2.50W1E-12, 2.50W1E- 
13 and 2.50W5 E-14) a s wel l as to the model with turbu- 
lence from Ri cher et al.l ( l2000l) which best reproduced the data 
(2.50T5.2D1M-4). Of the 16 observed elements which are in- 
cluded in our calculations, 12 (He, Li, O, Na, Mg, Si, P, Ca, 
Ti, Cr, Fe and Ni) are well reproduced by both the turbulence 
model and t he 2.50W1E-13 model. It is interesting to note that 
iBertin et alJ (fl995) determined, from Mgll lines, that the ob- 
served mass loss rate of SiriusA is between 1.5 x 10 and 
2x 10"'^. The overall fit is also just as good for the 2.50W5E-14 
model since it is also able to reproduce the carbon abundance, 
though it perhaps overevaluates the iron peak abundances. The S 
abundance is not at all reproduced by our calculatio ns; how ever, 
the observer gives little credibility to its value (HilL ll995h . It is 
clear that the model with a mass loss rate of 10"'^ Moyr^Moes 
not lead to the observed surface abundance pattern. Finally, sur- 
face abundance observations are not sufficiently accurate to en- 
able a differentiation betwee n turbulence and mass loss models. 

In iMichaud et al.l (l2005h . the binary system oLeonis (HD 
83808/83809) has been interpreted as consisting of two AmFm 
stars with masses of 2.12 and 1.87 Mq (Griffin, 2002). The au- 
thors show that two models with turbulent mixing, of 2.24 and 
1 .97 Mq respectively, are able to reproduce the observed features 
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Fig. 21. A model of 2.2 Mq with a mass loss rate of 5 x 
lO^'^^Moyr"' (solid line) and a model of 1.9 Mq with the same 
mass loss rate (dashed line) are plotted in an H-R diagram (a). 
The observed position in the H-R diagram of the primary and 
secondary components of o Leonis (with error bars) are shown 
with crosses. The squares indicate two possible positions on the 
primary's evolution path which fall within observational error 
bars. The circle indicates the position of the secondary at the age 
of the primary indicated by the square. The surface values of Fe 
and Ca are shown relative to their original values in (b) and (c). 



(namely both positions in the H-R diagram as well as surface 
abundances) of the A and B components of the binary system 
within the observational error bars. Similarly, models of 2.20 
and 1.90 with a mass loss rate of 5 x lO^'^^Moyr"' are able 
to reproduce the H-R position of both components. At the age 
indicated by the squares (750 Myr and 750.5 Myr), the A com- 
ponent is in the rapid evolution stage that follows the depletion 
of hydrogen in the core (the luminosity change is related to the 
star's adjustement to H-shell burning). At this age, the B com- 
ponent is still in the slowly evolving MS stage, which explains 
why both circles overlap on the graph. 

The two squares in Fig. |2T] show that the Fe abundance can 
vary between -0.2 dex and 0.6 dex within the T^s error bar for 
the A component, while the corresponding Ca surface abundance 
varies from -0.18 to -0.4. Throughout most of this interval, the 
2.20 Mo model has anomalies which are typical of Am stars (Fe 
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overabundance coupled with a Ca underabundance). For either 
of the values for the A component, the 1.90 Mq model has an 
overabundance of Fe of 0.4 dex coupled with an underabundance 
of Ca which attains -0.3 dex, both typical Am star anomalies as 
well. 

Both the turbulent model of iMichaud et al.l (l2005h and the 
mass loss model can reproduce the AmFm character of compo- 
nents A and B. While the AmFm character of component A can 
be fitted by the turbulent model for its exact observed T^s, the 
mass loss model only generates typical AmFm iron overabun- 
dances for a part of the error bar on the hot side of the observed 



7.2. Open Cluster Stars 

In Fig.|22] we compare our results to observed abundance 
determinations for the hot Am star 68Tau (Feff - 9050 K, 
iNetopil et al.l 120081) from the Hyades open cluster The cluster 
age has been q uoted between 625 Myr (jP erryman et al. 199^ 
and 783 Myr bv lVarenne & Monie^ (ll99S ). Its quoted metallic- 
ity has a lso ra nge d between Zq ^ 0.024 (Perrvrnan et alj.l998t 
lGrattonl l2000) Q and Zq = 0.03 (ICavrel et al.lll985h using F 
and G star iron abundances as indicators. In order to reflect this 
metallicity, the selected models have been calculated with an 
initial m etalhcity of Zp - .03, which was also used in mod- 
els from iRicher et all (|2000 H We have attempted to make a 
compromise between fitting age and Teff: three models have a 
mass of 2.50 Mq and one has a mass of 2.30 Mq. The 2.50W1E- 
13Z0.03 and 2.50W5E-14Z0.03 models are on a short T^ff up- 
swing which arises as hydrogen nears depletion in the core (in 
Fig.[2T|for instance, it is the segment which immediately fol- 
lows the main-sequence, spanning from logTeff - 3.87 at its 
bottom to logTeff - 3.93 at its top). In terms of stellar age, 
this upswing only lasts 3 Myr before the star starts its descent 
onto the red giant branch. While it has very little effect on sur- 
face abundances, models were chosen at this age in order to be 
closer to the star's surface temperature. Given the large spread in 
abundances between observers, the fit is almost perfect with the 
2.30W1E-13Z0.03 model, which is slightly cooler and younger, 
yet is still on the main-sequence. Of 15 observed elements, only 
Na, Al and Mn (arguably just Al) are not reproduced. The fit is 
as good if not better (because of Ni) than the fit obtained with the 
model with turbulence. The 2.50Wle-13Z0.02 model was added 
in order to illustrate the effect of reducing initial metallicity on 
absolute abundances (see also Fig. [13]). The fit with observations 
is better than for the 2.50W1-13Z0.03 model, since the iron peak 
abundances with respect to solar abundances are smaller in the 
lower metallicity model. By comparing these two curves one can 
conclude that, for a given mass loss rate, a 0. 18 dex reduction of 
initial metallicity can, at most, lead to a 0.1 8 dex reduction for 
elements which are not supported, such as C, N and O, and a 
0.09 dex reduction for elements which are supported such as Fe. 
See also Sect. l5331 

In Fig.|23] we compare 2 models of 1.9 Mq with mass loss 
as well as a model with turbulence to the observations of the 
star HD 73045 (Teff - 7500 K) from the Praesepe open cluster 



" Using the Hyades iron enrichm ent factor of iGrattonj (l200(]h to mul- 
tiply the metallicity determined by I Asplund et al.l ( 120091) foTthe solar 
mixture, the Hyades metallicity becomes Zq - 0.019. The actual value 
would likely lie somewhere between 0.019 and 0.03. The curve with 
Zo = 0.02 in Fig.|22]illustrates the impact of uncertainties. 

Of course, [N/H] at / = for models with Zq = 0.03 is above solar 
for all metals. 
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Fig. 22. Observed surface abundances of 68 Tau (also 
known as vB 56, HR 1389 or HD 27962), the hottest 
star (bl ue straggler) from the Hyades open cluster. Circles 



Hui-Bon-Hoa & Alecian _i 19981): triang les, Robv & LambertI 
(199^ 



squa res. 



Burkhart & Couory (Il989h : asterisks 
Takeda & Sadak ane (1997). communication. Calculated values 
are shown for 4 models with varying mass loss rates as well 
as the model with turbulen ce which bes t reproduced the data 
(2.30R1K-3Z0.03, see Riche r et al.ll2000h . One model was cal- 
culated with Zo - 0.02, while all other models were calculated 
with an initial metallicity of Zq = 0.03. Metallicity is indicated 
in the model name. 



which has an approximate age of 800 Myr and a solar metallic- 
ity. There are 15 observed elements which can be compared to 
our simulations, although 3 determinations (N, K and Mn), result 
from a single line and therefore could be inacurrate. Again, note 
the large discrepancies between observers. Only the 1.90W5E- 
14 and the turbulence model can reproduce either the overabun- 
dant iron peak elements or the underabundances of C and O. 
The abundances of Na and Si are not reproduced by either of the 
models. 

Finally, we have compared our models to observations of the 
Coma Berenices star HD 108486 (Fig.l24li. Coma Berenices is an 
open cluster with an age of about 500 Myr and with a metallicity 
which is about solar We have matched the star's T^ff and age 
quite well with two 1.8 Mq models with mass loss. Except for 
O and Na, which are not reproduced by any of the two models, 
most elements are fitted by both models. Assuming error bars 
for S and Al which are similar to those for other elements, we 
can state that 10 of the 13 abundances can be reproduced by the 
1.80W1E-13 model, and 9 by the 1.80W5E-14 model. 

Although we have opted not to add any extra figures, 
our results are also compatible with observations of Ca over- 
abundances (see Fig.fT2]l in very young open clusters such as 
the Pleiades (Gebran & Monier, 2008 . Hui- Bon-Hoa & AlecianL 
1998) and ffPersei (Hui-Bon-Hoa, 1999). This is noteworthy 
since models with tur bulence do not pre dict such an overabun- 
dance (see Figs. 10 of lRicher et ani2000h . 
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Fig. 23. Observed surface abundances of HD 73045 (Teff ^ 
7500 K) of the Praesepe ope n cluster which has an approxi - 
mate age of 800 Myr. Circles . iHui-Bon-Hoa & Alecianl (Il998h; 
square s, iBurkhart & Couprvl ( 20001) : iFossati et alJ 

(|2007|) . Curves correspond to models listed at the top of the 
figure. All models were computed with solar initial abundances 
(Zo = 0.02). 
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Fig. 24. Observed surface abundances of HD 108486 {T^s ^ 
8180K) of the Coma open clust er which has an ap- 
jroximate age o f 500 Myr. Circles, iHui-Bon-Hoa & Alecianl 



proxu ^ , _ . . 

(Il998i): squares. iBurkhart & Couprvl (l2000h: five-point stars, 
iGebran et al.l (l2008l) . Curves correspond to models listed at the 
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All models were computed with solar initial 
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Fig. 25. Lithium, beryllium, iron and calcium abundances for 
models with and without mass loss at 625 Myr, the approximate 
age of the Hyades open cluster All models were calculated with 
an initial metallicity of Zq = 0.02 and the original Li abun- 
dance was set at N(Li)=3.05. The Li observations are from (A, 
and 7 for upper- limits) Boesgaard & Tripicco (1986), ( a and 
T) IBoesgaard & Budge ( 1988)', (□) Burkhart & Couprvl (Il989h 
and (■) Burkhart & Coupry (2000). Be abundances are taken 
from (•) Boesgaard & King (2002). Additional calcium and iron 
abundances are also shown for all stars with lithium determi- 
nations: (x) Boesgaard & Friel 1990; (*) Takeda & Sadakaii^ 
ri997l: (o) iHui-Bon-Hoa & Alecian|[l998l . All stars with multi- 
ple determinations are connected by a line segments. Calculated 
models are indicated by dots along the curves. 



7.2.1. Lithium gap 

In Fig.|25] models with and without mass loss are compared to 
lithium, beryllium, calcium and iron observations in and around 
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the Hyades lithium gap. Lithium determinations are shown for 
F stars and AmFm stars (normal A stars and other peculiar A 
stars are omitted). Iron and calcium abundances are shown only 
for stars which had a li thium determination . All b eryllium abun- 
dances for F stars in iBoesgaard & KingI (|2002|) are showrQ 
When multiple observations for the same star were available, the 
different determinations are connected by a line segment. This 
gives an evaluation of the uncertainty. All models calculated 
with mass loss rates of 5 x lO^'^^Moyr"' and 1 x lO^'-^Moyr"' 
which were still on the main-sequence at 625 Myr are shown. 
Models with a mass loss rate of 1 x lO^'^Moyr"' are omitted 
since they result in surface abundances which are very simi- 
lar to diffusion only models (compare the diffusion only and 
1.40W1E-14 models in Fig.fTSTl. Therefore in the following dis- 
cussion, results from diffusion only models can be assimilated 
to models with unseparated mass loss < 1 X lO^'^^Moyr-'. All 
models were calculated with an initial metallicity of Zq = 0.02, 
although the metallicity of the Hyades is above solar (Z - 
0.024, see Sect. l7.2b . The original value of lithium was set to 
A(Li)=3.05. This value fits lithium determinations for the stars 
at the top of the cold side of the gap, which, since diffusion 
plays only a small role for these stars, probably reflect the clus- 
ter's original Li content (unless there is significant pre-main- 
sequence burning). Following the same logic, A(Be) was set to 
L40. 

According to our calculations, atomic diffusion in the ab- 
sence of competing processes leads to an important reduction of 
surface lithium abundance. The smallest lithium anomaly was 
obtained for the LIO Mq model, for which surface lithium was 
reduced by 0.015 dex at 625 Myr. The largest lithium reduction 
at 625 Myr is by about -L4dex (or a factor of 25), for the 
1 .46 Mq model witho ut mass loss. This does not qu ite reach the 
bottom of the gap of iBoesgaard & Tripiccol (ll986H^'l As seen 
in Fig. [16] the 1.47 Mq model without mass loss would likely 
have reached a lower lithium abundance had it been able to con- 
verge up to the age of the Hyades, thus reconciling some of the 
difference. In fact, the mass of the model which would have at- 
tained the lithium gap minimum can be deduced from Fig.|3] 
Because Li is not supported until gi-ad - g just below the sur- 
face convection zone at 625 Myr, the heaviest model for which 
the BSCZ is located where log T > 5.4 throughout its evolu- 
tion will represent the gap minimum, since it is for this model 
that Li is sinking fastest. From Fig.|3] while the 1.55 Mq model 
is clearly on the hot side and the 1.43 Mq model on the cold 
side, it is the 1 .46 M© model that should be closest to the bot- 
tom of the gap. Furthermore, if the competition with Fe that is 
prescribed by the OPAL opacities and calculated in these mod- 
els without mass loss is correct, then diffusion alone cannot ex- 
plain the increase of Li on the hot side of the gap. The grad(Li) 
in the atomic diffusion only models will remain smaller than g 
by a factor of at least 2. However, given the uncertainties dis- 
cussed in Sect. l3.1.T] on the location of Fe lines, grad(Li) could 
very well attain g near log T ^ 5.3, in which case lithium would 



L ithium abundances from iBoesgaard & Tripiccol il986h and 
Boesgaard & Budge ( 198^) are prefered over the revised values from 
Boesgaard & Kina 62002') because there is a greater number of stars. 
Nonetheless, in the more recent paper, all previous determinations are 
revised upwards by 0.09 to 0.4 dex. 
-° It has been suggested to revise these ob servations up wards by up 
to 0.4 dex (see the discussion in Sect. 2.1.2 of lMichaud & C harbonneau 
11991! and IBoesgaard & King|l2boa) . The maximum depletion encoun- 
tered in the deepest part of the gap might so be closer to a factor of 
50. 



be supported, and would consequently exhibit different surface 
behavior. 

For the models without mass loss on the cold side of the 
gap, the Fe abundances are in agreement with observations up 
to about 6800 K, after which the calculated overabundances be- 
come too large. The discrepancy between the Fe curve and the 
observations for Teff < 6700 K is related to our models having 
a solar initial metallicity, whereas the Hyades stars were formed 
in a metal rich environment. 

The models with mass loss rates of 5 x lO^'^^Moyr"' and 
10 '^Moyr ' cannot explain the depletion encountered within 
the gap. However, given reasonable error bars, they are con- 
sistent with the almost constant lithium and beryllium abun- 
dances observed for T^ff > 7200 K. Both mass loss rates lead to 
models which reproduce the observed Fe abundances between 
6000 < Tgff < 8200 K and, in particular, the increase in Fe abun- 
dance for Teff > 7200 K, which is compatible with the AmFm 
character of these stars. Given the large discrepencies in determi- 
nations, most Ca abundances are also compatible with the mod- 
els with mass loss. 

Neither the calculated gap minimum nor the shoulder on the 
cold side of the gap match the observed positio n in Teff . By com- 
parison to Fig. 9 of [Richer & Michau dl ([T993l) . their calculated 
depth for the gap (- 1 .6 dex) resembles the depth obtained in our 
calculations (-1.4 dex, see discussion in Sect. l5.3.4] l. The shoul- 
der on the cold side of the gap obtained in the present calcu- 
lations matches the curve they obtained with Zq = 0.02 within 
+50 K. Accordingly, as also seen in this same Fig. 9, if we had 
chosen Zo = 0.03 ([Fe/H]=-i-0.18), some of the 150 K difference 
would have been recuperated as the gap minimum would have 
been shifted toward cooler temperatures by 50-80 K. There is 
also a +50- 1 00 K uncertainty on the observed potition of the gap 
(see discussion in Sect. 2.1.2 of lMichaud & Charbonneaull99ll) . 
The uncertainty on the age of the Hyades (from 625 to 783 Myr, 
see Sect. 17. 21 ) could also account for some of the difference as 
illustrated in Fig.[T71 As the age of the isochrones increases, the 
Tss at which lithium abundances fall off also decreases. The real 
problem in explaining the Li gap with atomic diffusion is not 
with the exact T^ff of the gap nor its depth, but rather with the 
calculated Fe overabundances which are not observed, and the 
related difficulty in calculating §rad(Li) on the hot side of the gap 
due to Fe lines. 

8. General discussion and conclusion 

8.1. Summary of results 

Evolutionary models including both atomic diffusion and un- 
separated mass loss explain the main abundance anomalies of 
AmFm stars (Sect. l7.2] l. When mass loss is assumed to be the 
only macroscopic process competing with atomic diffusion, ob- 
served abundance anomalies from open cluster stars as well as 
Sirius A and o Leonis constrain mass loss rates to 2 - 5 times 
the solar mass loss rate. As shown in Sects. 7.1 and 7.2, models 
involving mass loss are as capable as models involving turbu- 
lence in explaining observations of AmFm stars. This is because 
in both instances, the important separation occurs at the same 
depth (AM/M» ^ 10"* - 10"^) for most of the main-sequence 
life. Whether the mass loss model is to be preferred over the tur- 
bulence model is difficult to assess given the large observational 
uncertainties. However, as shown in Fig. [191 the internal distri- 
bution of elements is different between the two cases for most 
elements. With differences reaching a factor of 4 - 5 for abun- 
dant elements such as Fe, there should be effects on local opac- 
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ities and thus on pulsations. Asteroseismic tests could perhaps 
distinguish between the twcPl 

In the mass loss regime, chemical separation affects up to 
10^^ Mq of a star's mass or, equivalently, 20 to 25% of the 
stellar radius (see Fig.|6] and Fig.[T9]l. For any given element, 
as long as the wind velocity is greater in amplitude than the 
downward settling velocity, the local abundance solution is de- 
termined by flux conservation; local abundances adjust as the 
flux quickly becomes constant throughout the outer envelope 
(see Fig.|7]i. As a result, the surface abundances depend on mat- 
ter which is advected from deep inside the st ar (see discussio n in 
Sect. l5T2] i. This diff^ers from the models of lWatsonI (1 197 l b and 
lAlecian {199^ in which surface abundances depend on the outer 
10"'° Mq. This can also be contrasted to the solution obtained in 
the models with weak or fully separated mass loss presented in 
Sects.EUJand 15331 

When flux is conserved throughout the envelope, abundance 
gradients which form near the surface, between surface convec- 
tions zones for example, have no effect on the surface solution 
once the star has arrived on the main-sequence (see discussion 
in Sects. l57J!2l and l5?3.3l l. In this instance, if the age of interest 
is greater than Mbscz/M, where Mbscz is the mass above the 
bottom of the deepest surface convection zone, one can obtain 
a nearly similar surface solution by approximating that abun- 
dances are homogeneous from the surface to the bottom of the 
deepest SCZ. However, early in the evolution, only matter from 
superficial layers has had time to be advected, and thus surface 
abundances obtained here depend on separation that occured 
close to the surface as first studied for Ca by Alecian! (1 19961 
see also Sect. l5.3.T] ) and confirmed observationally (see end of 
Sect. 17. 2b . This favors models involving mass loss rather than 
turbulence. Likewise, variations obtained near the surface, which 
do not appear in models with turbulence, have an effect on the 
PMS (Fig.O and will be discussed in a forthcoming paper. 

In all models heavier than 1 .47 Mq without mass loss or 
with an unseparated mass loss rate < 1 x lO"''*M0yr"', the ac- 
cumulation of Fe and Ni around T - 200 000 K leads to the 
appearance of a thin radiative layer which separates the iron 
peak convection zone from the surface H-He convection zone. 
This accumulation forms before the appearance of a small in- 
version of the local molecular weight gradient inversion (see 
Sect. l5.2.r|l. The inclusi on of thermohaline convection as sug- 
gested bv Thead o et al.l (l2009l) could have an effect on abun- 
dances in the region, though convection occurs even when there 
is no molecular weight gradient inversion. This would require 
further investigation. Nonetheless, the appearance (or not) of the 
iron peak convection zone does not have a significant effect on 
the surface solution, nor does it appear in models with mass loss 
which adequately reproduce observed abundance anomalies of 
AmFm stars (Sect.©. 

Since this paper is a part of a series which explores the vari- 
ous macroscopic processes which compete with atomic diffusion 
in AmFm stars, it is important that the models be as similar as 
possible to those used in previous calculations (e.g. those with 
turbulence) in order to isolate the effects due specifically to mass 
loss. This is one of the primary motivations for using the same 
initial solar abundances as in previous calculations, rather than 
the revised Asplund et al. (2005, 2009) abundances (see also the 
discussion in Sect|2]i. Varying too many things at once could 
obscure results and introduce further uncertainty. Furthermore, 
there is a controversy on solar abundances, since heliosismology 



strongly favors the older ( Grevesse et al.L Il996h over the newer 
(Asplund et al.L I2005L 12009 ) composition. One may then view 
the abundance differences between the two sets as an evaluation 
of uncertainty. Since solar abundances are used throughout this 
paper, the uncertainty on solar abundances leads to uncertainties 
on the absolute values of all abundances. As shown in Fig.l22l 
a factor of 1.5 reduction (or O.lSdex) of the original Z leads 
to a similar reduction (O.lSdex) of the expected abundances of 
atomic species that are not supported, as well as a smaller re- 
duction of 0.09 dex for species such as Fe which are supported 
by gr-dd- The fit for the abundances of 68 Tau is about the same 
for both values of Z as seen in Fig.|22] Equivalently, compen- 
sating the change of Fe abundance would require reducing the 
mass loss rate from lO^'-'Moyr"' to 7 x 10^'^Moyr"' according 
to the results shown on Fig.[TT] This may be viewed as the un- 
certainty on the mass loss rate resulting from the uncertainty of 
solar abundance^ 

8.2. Further implications 

Atomic diffusion alone cannot explain all characteristics of the 
Hyades lithium gap, nor can unseparated mass loss. The cold 
side of the gap can only be reproduced by diffusion only mod- 
els or models with M < 1 x 10 '"'Moyr"', whereas the hot side 
of the gap and the AmFm character of stars for Tetj > 7200 K 
require a stronger mass loss rate. Moreover, separated mass loss 
(see Sects. 15331 and 17.2.1b seems required to explain observed 
Li underabundances near the bottom of the gap as well as re- 
duce the calculated Fe overabundances. In Fig. [18] in compari- 
son to the diffusion only model, the curve for case 2 shows both 
larger underabundances of Li as well as smaller overabundances 
of Fe. Well tuned fully separated mass loss (case 1) could do 
the same. Similarly, the model for case 3 has a nearly flat Fe 
surface abundance coupled with similar Li underabundances to 
the diffusion only model. It does not seem justified to further 
speculate on the role of separated winds in Li gap stars until 
we have a better understanding of separation mecanisms within 
stellar winds. Since radiative forces generally increase with Tgff, 
the above mentioned increase in mass loss rate seems possible if 
winds of A and hot F stars are radiative in nature. Since a star's 
Teff changes over time, a mass los s rate w hich depends on Tpff 
(or on L,) could also vary in time dSwenso n & Faulkner, 199^. 
Such effects were not introduced in order to limit the number of 
adjustable parameters. 

The competition between atomic diffusion and meridional 
circulation in 2-D should lead to solutions which resemble 
those obtained with mass loss, since meridional circulation 
leads to an additional advective term in the transport equation 
(Eq.|7]i. Therefore, because the internal distribution of elements 
in the mass loss regime differs considerably from the varia- 
tions encountered in the turbulent mixing regime (see discussion 
in Sect.|6]l, internal distributions due to meridional circulation 
could also differ significantly from those encountered via turbu- 
lence. Hence, when building stellar models, one should be cau- 
tious when replacing meridional circulation, which is an advec- 
tive process, by turbulent mixing. A careful study of the atomic 
diffusion of metals within the context of meridional circulation. 



^' ICarrier et alj i2007h did not detect pulsations which could have 
been a signature of iron accumulation in the Am star HD209625. 



^- To significantly improve the evaluation of the effect of changing 
to another set of solar abundances would require first recalibrating the 
mixing lenght using a solar model, then carrying out calculations for 
AmFm stars for both turbulence and mass loss, as well as reanalyzing 
observations of AmFm star abundances using the new solar abundances. 
This is outside the scope of the present paper. 



26 



M. Vick, G. Michaud, J. Richer, and O. Richard: AmFm and lithium gap stars 



such as the one carried out for helium in lOuievv et alJ (l2009h . 
could help determine the implications of such an approxima- 
tion. This could perhaps lead to asteroseismic tests which could 
distinguish b etween m odels using rotationally induced turbu- 
lence (Talon et al.1'2006) and those using meridional circulation 
jCh arbonneau & Michaud 1988), which are both used to explain 
the disappearance of the AmFm character for rotation velocities 
greater than 1 00 km s ' . 

Observati ons of rapid p-mode oscillations in many Ap stars 
(iKurtz' T978I) an d in particular in Przybylski's star (see also 
[Mkrtichian^in 12008) have led to a number of studies of the 
oscillation mechanisms. In particular. IVauclair et al.l (1199 1^ sug- 
gested that unseparated mass loss acting solely in polar regions, 
where the magnetic field is strongest, could induce helium gra- 
dients which are compatible with oscillation generating models 
dBalmforth et al. 2001). A lthough differences be tween our stel- 
lar model and the one of lVauclair etal.1 (11991 ^ could have an 
effect on the predicted anomalies (notably the absence of con- 
vection due to magnetic braking/freezing in the latter), our cal- 
culations suggest anomalies would reach inwards to about 25 % 
of the star's radius and could have an important effect on opaci- 
ties. Though it would depend on the strength of the overall mass 
loss rate, which will be smaller than the mass loss rate at the 
poles, similar He depletions could be coupled with overabun- 
dances of iron peak elements around 200 000 K and perhaps iron 
convection. It is not clear whether magnetic braking/freezing 
or thermohaline convection can stabilize iron peak convection. 
Unfortunately, we are not able to investigate this scenario any 
further since our models require spherical symmetry. 

Perhaps asteroseismology will allow us to answer some 
of these questions, while revealing the relative importance of 
meridional circulation, turbulence and mass loss within chem- 
ically peculiar stars. 
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